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FINITE-DIFFERENCE THEORY FOR SOUND PROPAGATION IN A LINED DUCT 
WITH UNIFORM FLOW USING THE WAVE ENVELOPE CONCEPT 
by Kenneth J. Baumeister 
Lewis Research Center 

SUMMARY 

The finite-difference equations are derived for sound propagation in a two- 
dimensional, straight, soft-wall duct with a uniform flow by using the wave envelope 
concept. This concept involves a transformation of the wave equation into a form whose 
solution does not oscillate in the axial direction. This transformation reduces the re- 
quired number of finite-difference grid points by one to two orders of magnitude. In 
this report, the governing acoustic-difference equations in complex notation are derived. 
Also, a new exit condition is developed that allows a duct of finite length to simulate the 
wave propagation in an infinitely long duct. Sample calculations are presented for a 
plane wave incident upon the acoustic liner. They show the numerical theory to be in 
good agreement with closed-form analytical theory. Also, complete pressure and ve- 
locity printouts are given to some sample problems and can be used to debug and check 
future computer programs. 


INTRODUCTION 

Finite-difference techniques are now being considered for application in complex 
situations involving sound propagation in ducts with axial variations in Mach number and 
impedance, as well as variations in duct geometry. Before these more complex prob- 
lems are attempted, the finite-difference theory is applied to the simpler problems of 
noise propagation in ducts without flow and with steady uniform flow. 

In the absence of flow, numerical theories were successfully applied to straight 
ducts (refs. 1 to 4), variable- area ducts (ref. 5), and ducts with variations in the wall 
impedance (refs. 6 and 7). Most recently, numerical techniques were applied to noise 
propagation in turbofan variable- area inlets with hard walls (ref. 8). Reference 3 de- 
veloped the difference equations for noise propagation in an acoustically lined duct with 
uniform mean flow. Results were presented for a soft-wall duct in the absence of graz- 
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ing flow and for a hard-wall duct with grazing flow. However, reference 3 did not pre- 
sent any results for a soft-wall duct with mean flow because the soft-wall difference 
equations had not been completely programmed and debugged when reference 3 was pre- 
sented. 

For long ducts with relatively large attenuations, the theory and equations developed 
in reference 3 to simulate an infinite duct appear to be valid, based on a comparison of 
the finite-difference theory of reference 3 with the closed-form analytical solution of 
reference 9 (table xn). However, for short ducts with a length-height ratio less than 1, 
the numerical and analytical solutions were not in agreement. The transverse acoustic 
velocity calculated from the wall impedance condition was not equal to the same velocity 
calculated from the momentum equation. A large difference existed even though the two 
approaches should have yielded the same velocity at the wall. A number of possible 
causes were investigated, with the result that this problem finally was traced to the exit 
difference equation presented in reference 3. This report seeks to resolve the problem 
by deriving the exit condition that is valid for all duct length-height ratios. 

For long ducts or high-frequency sound propagation, many grid points are necessary 
for a proper numerical simulation. The large number of grid points necessitates large 
computer storage and long execution time. To determine the optimum liner impedance, 
many duct calculations are usually required in a search for the impedance that gives the 
highest duct attenuation. This considerably extends the computer running time. A wave 
envelope concept was developed (ref. 6) that greatly reduced the computer running time 
and storage requirements. 

In this report the wave envelope concept is applied to the case of uniform flow in 
lined ducts. The wave envelope concept was first developed in references 4 and 6 for 
noise propagation in ducts in the absence of flow. Kaiser and Nayfeh (ref. 10) later ex- 
tended this technique to include wave propagation in variable- area ducts, again in the 
absence of flow. 

The wave envelope concept attempts to remove most (if not all) of the axially oscil- 
latory part of the wave pressure profile and thereby to reduce the number of grid points 
required for a finite- difference solution. Rather than determining the actual pressure, 
the wave envelope concept transforms the governing equations such that they describe 
the envelope of the pressure. The number of grid points can be reduced by one order of 
magnitude. Reduction of grid points by two orders of magnitude is possible for problems 
with large duct length-height ratios and high frequencies. 

The derivation for the complete set of difference equations in complex notation is 
presented, as well as many of the details in the matrix solution. The complex form of 
the difference equations is simpler to program than the real form used in reference 3. 
These details are included so that the numerical techniques can be conveniently pro- 
grammed. Sample calculations are presented, and complete pressure and velocity print- 
outs are given to some sample problems that can be used to debiig and check computer 
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prc^rams. The limitations of the theory are also discussed. 


GOVERNING EQUATION AND BOUNDARY CONDITION 

The propagation of sound in a duct is described by a solution of the continuity and 
momentiun equations with the appropriate impedance boundary conditions. In this sec- 
tion, the basic governing equations and boundary conditions for a two-dimensional rec- 
tangular duct are presented. In the next section, the difference form of these equations 
is developed. 


Continuity and Momentum Equations 

The linearized equations for mass and momentum conservation can be written for 
a Cartesian coordinate system in the following form: 

Continuity: 


at* ax* ^ ^ \ax* ay*/ 


= 0 


Axial (x) momentum: 


au* = . _L 

at* ax* pq ax* 


Transverse (y) momentum; 


( 1 ) 


( 2 ) 


+ u* ^ (3 

at* ax* ay* 

where an asterisk denotes dimensional quantities. The assumptions involved in the der- 
ivations of equations (1) to (3) are given in reference 11, for example. (All symbols 
are defined in appendix A. ) 

As customary for steady state, the solutions for the dependent pressures and ve- 
locities are assumed to be of the form 


P*(x*, y*, t*) = P’(x*,y*)e^'^*** 


(4) 
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u*(x*, y*, t*) = U’(x*, y*)e^^*** 

(5) 

v*(x*, y*, t*) = v’(x*, y*)e^‘^*^* 

(6) 


where a prime denotes dimensional quantities with the time-dependency removed). The 
final solution will be represented by the real parts of the preceding quantities. Substi- 
tutii^ equations (4) to (6) into equations (1) to (3) and introducing the dimensionless pa- 
rameters 


P = 


(7) 


u = 




-P* 


M = 


^0 


( 8 ) 


X* y* 

X = — ; y = ^ 

H* H* 


(9) 


yield the dimensionless steady-state conservation equations 


P - 


i 

2nt]l 9x 




{2nrir 


2 \9x 9y 


( 10 ) 


u = i 9u 

0x 2 ittj 9x 


( 11 ) 


V = i I IM 9v 

9y 2nr} 9x 


( 12 ) 


The dimensionless frequency tj is given as 


7 ] =■ 


CO* H* f*H* 


2n c 


$ 




( 13 ) 


The dimensionless height y raises between 0 and 1, and the dimensionless length x 
between 0 and the dimensionless duct length L*/H*. 
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Wave Equation 


Equations (10) to (12) yield the dimensionless wave equation for vmiform flow: 


(1 - uh 


3xV 


9y^ 


i47rr7M-^ = 0 
dx 


( 14 ) 


Equation (14) in difference form could be solved to determine the pressure in the duct. 

In the exponential notation displayed in equations (4) to (6), the dimensionless pres- 
sure and velocities in general have both real and imaginary parts. Thus, 


P(x,y) = P^^^(x,y) + iP^^\x,y) 

(15) 

u(x, y) = u^^)(x, y) + iu^^^(x, y) 

(16) 

v(x, y) = v^^\x, y) + iv^^\x, y) 

(17) 


Therefore, each term in the wave equation (eq. (14)) is complex and requires special 
consideration in obtaining a solution. 


Wave Envelope Concept 

To remove most (if not all) of the axially oscillatory part of the wave pressure pro- 
file, the pressure P is transformed into a new variable p that chaises slowly in the 
axial direction. This new pressure p is defined as 

P(x, y) = p(x, y)e<-I2'A*)x = p(x,y)el-' 2 '’>V(l+M)l>‘ (18) 

where P is the acoustic pressure and p the pressure of the wave envelope. The real 
parts of P and p are shown in figure 1. The effective axial wavelei^th parameter X"*" 
is the wavelei^th of the oscillation in the soft-wall duct represented by figure 1. The 
wave envelope frequency parameter is the effective frequency of this pressure wave 
and is related to X"^ as 


_ 1 + M 

+ 

n 


(19) 
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This relationship was assumed because in a hard-wall duct, X'*’, Mach number M, and 
tj'*’ are exactly related as in equation (19) for plane-wave propagation. 

bi a soft-wall duct with sound propagating at a frequency of t], the effective fre- 
quency which would remove the oscillation according to equation (18), is not known. 
However, work performed in references 4 and 6 indicates that can be approximated 
by rj. Therefore, in most problems, tj'*' will be chosen equal to tj. However, since 
this may not always be the case, 77 "^ is not set equal to tj; rather, it is carried along in 
general form. 

Substituting equation (18) into the wave equation (14) yields a new governing differ- 
ential equation caUed the generalized wave envelope equation: 


K, + K„p + K, -^ = 0 

5x2 3 y 2 2 ax 


( 20 ) 


where 


= 1 - 


Kg = (2 ;t77)" 


1 _ (1 - M)| 


2M rj'^ 

(1 + M)l 


1 + M 77 


^3 = 


-Mn-rj |m + •^ (1 - M)j 


( 21 ) 

( 22 ) 

(23) 


Examples are presented to illustrate how equation (20) can remove the oscillatory nature 
of the pressure and thereby reduce the number of grid points necessary for a finite- 
difference solution. 

In soft-wall ducts, X"'" is not known precisely; therefore, the problem of picking 
X"*" (or T]"*") to remove the pressure oscillation must be considered. The uniform hard- 
wall duct, for which the correct answer is known, will now be investigated to determine 
the sensitivity of the final answer to the assumed value (guess) of rj'^. 

For plane waves propagating in hard-wall infinite ducts in the absence of flow, the 
wave envelope equation reduces to 


+ (2jr)^(rj^ - - i4j7T)‘^ = 0 

3x2 V y ax 


(24) 
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The exact analytical solution of equation (24) for the wave envelope pressure p is 

p = p^^) + ip^^) = cos 2Tr{r] - tt^)x - i sin 277(77 “ (25) 

The product of equation (25) with ^ gives the exact pressure fluctuation as de- 

fined by equation (18): 

P = pe"^^”^^ ^ + iP^^^ = cos 27rr?x - i sin 2nr}x (26) 

which is independent of the choice of iff", by definition. 

The anal 37 tical solutions for p^^^ from equation (25) for various assumed values of 
T]"*" or 1 77 - 77 ’*' I are shown in figure 2. For = 0, equation (25) reduces directly to 
equation (26); that is, p^^^ equals P^^^, see the cosine- shape line in figure 2. For this 
assumption, about 12 grid points would be required to describe this pressure profile 
adequately in a difference analysis. However, if, for example, rf^ is assumed to have 
a value between 0. 7 and 1. 3 such that |t 7 - t 7 '*"| < 0. 3, the exact anal 57 tical solution for 
p(^^ would require as few as three grid points in the difference analysis. The curves 
for 77 ’*' between 1. 0 and 1. 3 are the mirror image of those between 0. 7 and 1. 0, and 
therefore are not shown. K rf^ is assumed to be equal to 77 , the curve for p^^^ is the 
straight line t/"*" = 1, or 177 - 77 '^! =0. This line represents the envelope of the pressure 
oscillation. If, for example, it is assumed that 7 }'*' = 0. 8 , the exact solution in fig- 
ure 2 does not represent the true envelope of the pressure oscillation; it is still a gently 
varying function that requires only a few finite-difference grid points to describe its 
shape accurately. Thus, it is necessary only to pick a value of X'*' (or 77 ''') in the vicin- 
ity of the true value of X"*" to get a savings in the grid points required for a finite- 
difference analysis. In so doing, the differential Helmholtz equation is transformed to 
a new form (eq. ( 20 )) that requires fewer grid points in its solution. 

As another example of how to pick 77 ''', consider the problem of predicting the at- 
tenuation in a soft-wall duct, where ri"^ is an unknown. Calculations for the attenua- 
tion were made in reference 4 by assuming that 77 "^ = 77 . Figure 3 (from ref. 4) shows a 
comparison of the analytical and wave envelope calculations for the optimum soft-wall 
duct attenuation for various 77 and L*/H* values. 

Excellent agreement between the analytical and numerical calculations was obtained 
for the two-dimensional duct example in reference 4. For the 77 = 5 and L*/H* = 6 
case in figure 3, the conventional finite-difference theory required 3600 grid points; the 
wave envelope difference theory required only 100 grid points. Thus, a savings of 
3500 grid points over the conventional difference theory was obtained. At lower 77 and 
L*/H*, the savings was smaller. 

In principle, any choice of 77 ’*' will yield the correct answer; however, a poor 
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choice of rf*’ will require more grid points than a good choice. Thus, for any assumed 
value of rf", it is necessary to check for a converged answer by increasing the number 
of grid points in both the x and y directions, as was done in reference 4. 


Wall Boundary Condition 

The boundary condition at the surface of a soft -wall duct requires that at the wall 
the acoustic displacement of a particle in the fluid just outside the soft wall be equal to 
the displacement of a particle just inside the soft wall (ref. 11). In addition, the pres- 
sure and velocity fields at the boundary can be related in terms of the specific acoustic 
impedance where 


jf = = 2nnP 

w 




(27) 


As a consequence of these two relationships, the fluid velocity at the wall is 


vl 

C V 2nv 3x / 


(28) 


w 


and the transverse pressure gradient at the wall is 


dP 


-i2}T7]P 

2M 8P 

^ iM^ d^P 

9y 

W 


? 0X 

w ? 2 iri] 9^2 


(29) 


The derivations for equations (28) and (29) are given in appendix B, equations (B12) 
and (B13), respectively. 

Substituting the expression for the transformed pressure (eq. (18)) into equation (29) 
3 delds 


9y 


w 




w. 



w 


(30) 


where 
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w„ = (33) 

?2;t77 


Exit Boundary Condition 

The exit impedance must be specified at either the duct exit or some point in the 
far field. The exit impedance presented here will allow the numerical solution to ap- 
proximate the analytical results for wave propagation in an infinitely long duct. The 
wave propagation in the entrance region of an infinite duct can be represented by a finite 
duct length L* by choosing the exit impedance at L* so that no reflections occur. For 
a plane wave propagating in a hard-wall duct with uniform mean flow, the condition for 
no reflections at the duct exit is an impedance of PqCq, or 


(34) 


In the limiting case of zero flow (refs. 2 and 4) and for uniform flow in long ducts, cal- 
culations indicate that equation (34) is a reasonable approximation for soft-wall ducts. 

The axial pressure gradient at the exit can be related to the pressure at the exit by 
using the exit impedance condition to eliminate the acoustic velocities from the boundary 
condition. This is necessary since the exit condition is to be used in conjunction with the 
wave equation, which contains only pressure terms. Combining equation (34) with the 
continuity equation and the two momentum equations (see appendix C for derivation) 
yields the following exit condition (eq. (C15)): 
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e 


( 35 ) 


9x 





m(i + M \ 

\ ^e/ 8^P 

2ff?7(l - M^) 0y^ 


e 


Reference 3 used the exit equation for an infinite, one -dimensional, hard-wall duct; 
consequently, the second term in equation (35) was neglected. As mentioned in the IN- 
TRODUCTION, the failure of the earlier formulation to predict accurate transverse 
velocities in short ducts, L*/H* < 1, was traced to the neglect of the second term in 
equation (35). 

Again, substituting the expression for the transformed pressure (eq. (18)) into equa- 
tion (35) yields 


where 


0X 


= Ljp 


ay' 



(36) 


(37) 


(38) 


In the sample calculations to follow, will be taken as 1. 

Because of the nature of the finite-difference approximations to be presented later, 
equation (30) caimot be used at the exit of the duct. At the exit corner of the duct, the 
transverse gradient at the wall is given as (eq. (D39), appendix D) 


0P 


9y 


Iwall exit 



M 0P 



wall exit 


(39) 
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Equation (39) is a new equation; it was not previously developed. Equation (39) includes 
both the wall impedance ^ and the impedance at the exit The derivation of equa- 

tion (39) and its full rationale are discussed in detail in appendix D. 

Substituting the expression for the transformed pressure (eq. (18)) into equa- 
tion (39) yields 

= ^Njp + Ng (40) 

wall exit ^^'^wall exit 

where 



Entrance Conditions 

At the entrance, a uniform pressure profile 

P(0, y) - 1 (43) 

is used. A more general entrance condition, where P(0, y) depends on y, is equally 
easy to treat. In terms of the envelope pressure p, equation (43) becomes 

p(o, y) = 1 (44) 

Axial Acoustic Power 

The sound power that leaves a duct and reaches the far field is related to the in- 
stantaneous axial intensity at the duct exit. On the basis of the discussion in references 
12 and 13, the intensity can be expressed as 
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mil I I II I 


1 1 II III II II 


^ = (1 1 m2)P%* + + Mc5p5>i*2 


(45) 


This equation is an improvement on Ryshov and Shefter*s intensity equation (ref. 14), 
which was used in reference 3. The time-averaged intensity over a cycle is 


1^ = — 1^ = Real f — ^ Pu^ + MPP + — M_uu 


Pt' 


2irr] 


(46) 


{2irvr 


where the bar denotes the complex conjugate. 

Assuming that the axial velocity can also be expressed in wave envelope form as 


u u eL-^2Trrf' /{1+Mj]x 


(47) 


and using equations (18) and (47) give the expression for the time-averaged intensity as 


= Real|^-i^ pu^ j + Mpp + — — u^u^^ 


2717] 


(48) 




Thus, the formula for intensity in terms of p and u„ has the same functional form as 

ci 

equation (46). 

The total dimensionless acoustic power is the integral of the intensity across the 
test section 




(49) 


By definition, the sound attenuation (the decrease in decibels of the acoustic power from 
X = 0 to x) can be written as 


A dB = 10 Ic^iol 


/E 


Er 


(50) 
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Pressure Amplitude and Soxmd Pressure Level 


The pressure amplitude | Pj is a useful quantity to indicate how the acoustic pres- 
sure varies in the duct. By definition, 


P 



In terms of the wave envelope pressure p. 


P 



(51) 


(52) 


FINITE-DIFFERENCE FORMULATION 

Instead of a continuous solution for pressure, the pressure is determined at iso- 
lated grid points by means of the finite-difference approximations, as shown in figure 4. 
Thus, the differential equations can be changed to a system of algebraic equations in 
terms of the pressure at each point. This set of algebraic equations can be solved 
simultaneously to determine the pressure at each grid point. 


Difference Equations 

The governing difference equations can be developed by an integration process in 
which the wave envelope equation (eq. (20)) is integrated over the area of the cells shown 
in figure 4. 


j j -^-2 + -^-2 + KgP + Kg jdx dy = 0 
Cell area ^ ' 


(53) 


The finite-difference approximations for the various cells shown in figure 4 are ex- 
pressed in terms of the coefficients 

Vi, j+i + Vt+1, j “ “ 

The subscript k denotes the cell number. The derivation and expressions for the co- 
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efficients are given in appendix D. These coefficients are listed in table I. 


Matrix Solution 

The collection of the various difference equations at each grid point forms a set of 
simultaneous equations that can be expressed as 


IaICe] = [F] 


(55) 


where {A} is the known coefficient matrix, [p] is the unknown pressure vector, and [f] 
is the known column vector containir^ the various initial conditions. The matrix is 
complex. 

To illustrate the detailed structure of the matrix given by equation (55), consider the 
example shown in figure 5. For this case the detailed matrix structure becomes 


Cg dg Ug 

1 1 1 ~ 


'Plf 



bi Cj dj j ej 

bi Cj dj 1 ej 

1 \ i \ \ 


P12 

Pl 3 


-1 

^2 ^2 1 

®2 1 j j 


Pl 4 


-^2 




P2I 



ai jbj Cj dj 

1 \ i \ 


P22 


0 

H 1 ^1 ^1 

“li '1 1 ^ j ^ 


P23 


0 

^2! ^2 

^= 2 1 I 1 


P24 


0 

1 ^3 

1 *^3 1 ®3 j 


P3I 



\ 1 ^1 
\ 1 

1 ®i j j \ 

1 bi Cl di ) 6 i 1 \ 


P32 

P33 

= 

0 

0 

1 

1 

^2 1 *^2 ®2 1 ®2 1 


P34 


0 

1 

1 

\H 1 ®3 ®3 

1 1 I 


P4I 


~ 0 ~~ 

\ \ 

H S ®i ^1 1 ®i 


P42 


0 

\ 1 \ 

1 ^1 I ®i ^1 1 ®i 


P43 


0 

1 

^2 1 ^^2 *^2 ^ ®2 


P44 


0 

1 

1 

1 1 ®6 ^6 


P5I 


0 

\ j \ 

1 \ 1 ^4 1 64 C4 d 4 

j \ j ^4 1 b4 C4 d4 


P52 

P53 


0 

0 


1 1 ^ 5 | ’’S ® 5 _ 


_P 54 _ 


0 


(56) 
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where the derivation of the coefficients is given in appendix D, The dashed lines illus- 
trate where the matrix can be partitioned into a form that is*tridiagonal, which in more 
abbreviated notation becomes 



This can be generalized for any number of points as 



( 58 ) 
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The difference equations are cast in this form and solved by using the computer pro- 
gram written by Quinn (ref. 5), which directly inverts the matrix. The solution of the 
matrix gives the pressure at each of the grid points. These pressures can then be used 
to determine the acoustic velocities and intensities throughout the duct. 


ACOUSTIC PARTICLE VELOCITIES 

The axial and transverse momentvun equations are solved to determine the acoustic 
velocities in the duct. In turn, these velocities, in conjunction with the acoustic pres- 
sure, are used to determine the intensity levels in the duct and the attenuation that the 
duct provides. 


Axial Velocity 

The axial momentum (eq. (11)) can be rewritten as 

— + i^u = -^ (59) 

0X 0X 

where 

^ (60) 
M 

In terms of the transformed pressures and axial velocities, substitutii^ equations (18) 
and (47) into equation (59) yields 

— 5 + ijS /l - = -/3 -^ + i/3crp (61) 

dx \^/ dx 


where 


1 + M 


(62) 


Runge-Kutta solution . - Since the axial pressure gradient in equation (6 1) can be ob- 
tained directly from the numerical solution, equation (61) was treated as a simple or- 
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dinaxy differential equation. As a result, a standard fourth-order Runge-Kutta integra- 
tion was used to solve equation (61). The initial condition used for the solution was 
found by rewriting equation (34) for = 1 as 


u. 


"e = 

= 2ffTjP I , 


(63) 

(64) 


The Runge-Kutta technique worked well when the wave envelope technique was not used; 
however, this technique became luistable in many cases when the wave envelope tech- 
nique was applied. The instabilities resulted from the large steps in the independent 
variable x that were used in conjunction with the wave envelope approach. To correct 
this problem, an exact solution to equation (61) was employed. 

Integral solution . - The problem associated with the Runge-Kutta solution could be 
circumvented by using an integration technique. Solving equation (61) and applying 
boundary condition (64) yield 


u_(x, y) = 27T77 P(x , y)e 


i«p(x--x) 


+ /3e 


■icpy. r 

•/x 


•e ^icps 3p(s, y) 
ds 


- i/Sue 


/ X i 
e gi 


p(s, y)ds (65) 


where 


(p =p - (J 


( 66 ) 


Equation (65) has a disadvantage in that the pressure gradient is required for the inte- 
gration process. However, equation (65) can be further simplified by integration by 
parts. Recognizing that 



[e^^® p(s, y)]j ds = e ® p(x^, y) - e^^^ p(x, y) 


(67) 


and that 
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ei<!^sMSiJ^ds = e 

9s 


i<px^ 


P(Xg, y) - p(x, y) - icp 


/■ 


e^'^®p(s,y)ds (68) 


and substituting equation (68) into equation (65) result in the following expression for 
the axial velocity: 


ujx, y) = 2ni]^ + -^e ® ^ p(Xg, y) - /3p(x, y) - i/3^e“^^^^ 


p(s, y)ds (69) 


The advantage of equation (69) over equation (65) is that the integrand can be found di- 
rectly from the pressure field. 

The velocity at any point can be found from equation (69) by using the values of the 
pressure field from the solution of equation (55). Evaluatii^ the integral in equation (69) 
can be a problem. The parameter <p in equation (69) can be written 


M 



1 + mU/ 


(70) 


Therefore, the integral term in equation (69) can be written as 



e^^® p(s, y)ds = 


-r 


p(s, y)exp 



1 + M TJ / 


ds 


(71) 


The exponential in equation (71) represents an oscillating function whose wavelei^h can 
be considerably shorter than the wavelength associated with the pressure field. In fact, 
in the limit of small Mach numbers, the wavelength associated with the exponential in 
equation (71) approaches zero, but the wavelength of the pressure wave is simply the 
inverse of rj. 

To illustrate this more vividly, consider the example of a soft-wall duct with a 
dimensionless frequency rj of 1, an L*/H* of 3, and a Mach number of 0. 5. A 
typical pressure profile in a suppressor duct is shown in figure 6(a) by the dashed 

lines; the solid lines denote the amplitude of the exponential e^^^. Figure 6(b) shows 
similar results when the wave envelope concept is used. In general, as well as in this 
example, the grid spacing in the finite-difference theory is chosen sufficiently small so 
that the pressure field can be accurately determined. The grid spacing associated with 
the accurate determination of the pressure field will usually be too coarse to adequately 
evaluate the integral in equation (71) by a conventional numerical integration formula 
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such as Simpson's rule. Using a larger mmiber of grid points would be wasteful, in 
general and, in particular, would defeat the basic idea of the wave envelope concept, 
which is to reduce the number of grid points. Therefore, a combined analytical and 
numerical technique was used to circumvent this problem. 

Polynomial- e xponential integration . - The integral in equation (69) was expressed 
as the sum of the integrals between each grid point to give 


f ® p(s, y)e ds = 

Jk 


i ^i 


p(s, y)e^'^® ds 


(72) 


A fourth-order polynomial fit of the pressure p was constructed from the solution of 
the wave equation for every value of i. In turn, equation (72) could be integrated 
exactly. The exact formulas are given in appendix E. This technique proved very suc- 
cessful. 


Transverse Velocity 


Once the axial velocity is known, the transverse velocity can be quickly found by 
using the condition of irrotationality 


9u _ dv 
9y 9x 

which is developed in appendix C, equation (CIO). 

Substituting equation (73) into equation (12) gives 

V _ i _9P ^ iM 9u 
9y 2nrj 9y 


(73) 


(74) 


Thus, the solution of the transverse velocity can be found by a simple differentiation 
process once u is known. As a check, the solution of the transverse velocity can also 
be found by solving the continuity equation (eq. (10)) and applying the boimdary condition 
at the wall given by equation (28). This check led to the discovery of the problem dis- 
cussed in the INTRODUCTION. 

In terms of the wave envelope pressures, equation (74) remains unchanged, as 
given by 
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(75) 


V 


a 


iiE + iM.^ 
3y 2nr] 9y 


AXIAL INTENSITY AND POWER 

In terms of the difference notation, the axial intensity given by equation (48) can be 
written as 


I. . = Real[i±^p. m ^ .) + Mp. ^ + — u j -u . . 

y 2irn hi hi 27rrf ^yhl 


(76) 


The total power across a particular cross section, as given by equation (49), is written 
in difference notation as 



(77) 


By evaluating E^ at the entrance and exit positions, for use in equation (50), the sound 
attenuation for the duct is determined. 


SAMPLE PROBLEMS 

As the first example of the difference technique, consider the problem discussed in 
the INTRODUCTION. Recall, for small L*/H*, the transverse velocities as calculated 
from the momentum equation and boundary condition were not equal, even thoi^h the two 
approaches should have yielded the same velocity at the wall. An example of this dis- 
crepancy is shown in figure 7(a). Now, however, for the new exit condition (eq. (35)), 
the transverse velocities calculated from the momentum equation (eq. (74)) and the 
boundary condition (eq. (28)) are in good agreement, as shown in figure 7(b) . This 
agreement comes about because the pressure field has changed due to the introduction 
of the transverse gradient in the exit condition (eq. (35)). The attenuation calculated by 
using the exit condition of reference 3 is 6. 1; the attenuation calculated by using equa- 
tion (35) as the exit condition is 3. 4, which is a significant difference. 

In another example of the difference technique, the noise attenuation at the optimum 
point (point of maximum attenuation in the impedance plane) is calculated for a two- 
dimensional duct with an L*/H* of 3. 43, a Mach number of 0. 3, and the ranges of the 
parameter 77 and ? listed in table n. In figure 8, the numerically calculated attenua- 
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tions are compared with corresponding analytical results (ref. 9) that are applicable to 
infinite ducts. The numerical results (lines) and the analsrtical results (circular ssrm- 
bols) for the maximum attenuation are shown in figure 8 to be in very good agreement. 

For symmetrical ducts (impedance identical at both walls), the pressure profiles 
are symmetrical about the centerline of the duct. In this case, 9P/9y can be set equal 
to zero at the centerline of the duct. Consequently, the number of grid points in the 
y-direction could be cut in half. 

The deviation between the closed form and the numerical analysis shown in figure 8 
at the lowest dimensionless frequency is assumed to occur because the PqCq exit im- 
pedance used in the ntnnerical analysis begins to deviate from the exit impedance in the 
anal 3 d:ical results (ref. 15, p. 246). By extending the soft-wall duct (ref. 2, appendix E), 
it is possible to reduce the reflected contribution and thereby more closely simulate an 
infinite duct, which is the basis for the analytical calculations. The dashed line in fig- 
ure 8 results from increasing the duct L*/H* to 6 . 86 and evaluating the attenuation at 
an X of 3. 43. 


PROGRAMMING AIDS 

Appendix F contains the solutions of four additional sample problems. Complete 
pressure and velocity printouts are given and can be used to debug and check the com- 
puter programs. The pressure profiles for the two standard solutions that do not use 
the wave envelope concept = 0 ) were checked by two independent computer programs 
(refs. 2 and 5) that used two different techniques for solving the matrix equations. 
Therefore, if the same difference equations are programmed, agreement of the pressure 
field to three-place accuracy should be expected. 

The outputs of the velocity fields are also presented in appendix F. For the 
uniform-velocity cases, these flow fields have not been double checked by independent 
programming because these solutions have been developed from the new techniques pre- 
sented in this report. 


LIMITATIONS OF THEORY 

At low frequencies (77 ^ 0. 5), the exit impedance PqCq used in this analysis does 
not necessarily lead to a condition of no reflection at the duct exit. Therefore, if the 
reflections are large, the numerical solutions will not necessarily be equal to the ana- 
Isdical infinite-duct results. For example, the attenuation calculated by the analytical 
infinite-duct theory for figure 7 was 2. 0 decibels, which is 1. 4 decibels lower than the 
numerical results. Doubling and tripling the duct length produced large oscillations in 
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the pressure amplitude |p| in the axial direction. This oscillation was interpreted to 
mean that large reflections were occurring at the exit of the duct. These calculations 
were performed without using the wave envelope theory. At low frequencies the wave 
envelope theory is not required since the wavelength of the pressure is very loi^. 

In general, a pressure wave is composed of a forward-goii^ and a reflected wave. 
The present wave envelope theory is useful only in those cases where the reflected wave 
is small compared with the forward-going wave, since the reflective component has been 
neglected in equation (18). The wave envelope technique may also be limited to the case 
with a single dominant mode or to the special multimodal case where the modes have 
similar effective axial wavelengths. In all cases, in solving a particular problem by 
using finite-difference theory, the grid spacing should be doubled in order to check for 
convergence. 

Finally, assuming equals 1 in equation (34) applies strictly to a single plane 
wave. Therefore, some reflection will occur at the duct exit for pressure waves other 
than plane. However, for tj > 1, this reflection is generally small. 

CONCLUDING REMARKS 

A finite-difference wave envelope theory for sound propagation in a two-dimensional 
soft-wall duct with uniform mean flow has been presented. Previous difference theory 
did not give valid attenuations or acoustic velocities for ducts with L*/H* of the order 
of 0. 5. The difference theory with plane-wave input developed herein is shown to be in 
good agreement with closed-form analysis for the complete rai^e of duct parameters. 

Also given are the latest numerical procedures in appl 5 dng finite-difference tech- 
niques to ducts. The derivations for the various procedures are presented in appen- 
dixes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 6, 1977, 

505-03. 
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APPENDIX A 


SYMBOLS 


A 

coefficient matrix 


submatrix 


matrix element 


submatrix 

\ 

matrix element 

C 

constant 

Ci 

submatrix 

c* 

speed of sound, m/sec 


matrix element 


matrix element 

AdB 

sound attenuation, dB 

E. 

1 

total acoustic power (eq. (77)) 

E 

X 

acoustic power (eq. (49)) 


E^ evaluated at x = 0 

*^k 

matrix element 

F 

initial- condition column vector 

^1 

initial- condition submatrix column vector 

f* 

frequency, Hz 

H* 

channel height, m 

Ix 

acoustic-intensity axial direction 

i 


Ki,K2,K3 

coefficients of generalized wave envelope equation (eq. (20)) 

L* 

length of duct, m 

Li, Lg 

exit- condition coefficients (eq. (36)) 

M 

Mach number, U*/cq 

m 

total number of grid rows 

N 

index: 0, 1, or 2 



wall exit coefficients (eq. (40)) 
total number of grid columns 

dimensionless pressure fluctuations P(x, y) = P*/P^ 
submatrix- colvunn pressure vector 
pressure fluctuation, P*(x*, y*, t), N/m^ 
amplitude of pressure fluctuation, N/m^ 

P’(x*, y*), N/m^ 

pressure- column vector (eq. (55)) 

wave envelope pressure 

element of pressure at i, j grid point 

polynomial coefficients 

dummy variable of integration 

time, sec 

mean flow velocity in x-direction, m/sec 
dimensionless acoustic particle velocity in x-direction 
dimensionless axial velocity (eq. (47)) 
dimensionless acoustic particle velocity in y-direction 
dimensionless transverse velocity (eq. (75)) 
boundary- condition coefficients (eq. (30)) 
dimensionless axial coordinate, x*/H* 
starting grid position for polynomial fit 
axial grid spacing 

dimensionless transverse coordinate, y*/H* 
transverse grid spacing 

O 

acoustic impedance, kg/m -sec 
dummy variable of integration 
coefficients of polynomial fit 
27r?7/M(eq. (60)) 

1/i^ Ax (eq. (E16)) 
specific acoustic impedance 



V 


dimensionless frequency 
wave envelope frequency 
effective axial wavelength 
particle displacement 

O 

density of uniform medium, kg/m 
2irr;V(l + M) (eq. (62)) 
potential function 
j3 - cr (eq. (66)) 
circular frequency 
Subscripts: 
e exit condition 

i axial index (fig, 4) 

j transverse index (fig. 4) 

k cell index 

opt optimum 

w wall 

X axial position 

y transverse position 

Superscripts: 

* dimensional quantity 

* dimensional quantity with time function removed 

( ) complex conjugate 

(1) real part 

(2) imaginary part 


+ 

V 

X+ 

^0 

Q 

# 

(p 

<JL>* 



APPENDIX B 


DISPLACEMENT CONDITIONS 

This appendix develops the relationships between the transverse velocity just inside 
an acoustic absorber and the transverse velocity just outside the absorber in the fluid 
medium. The development follows that of reference 11. By definition, the acoustic 
velocity of the medium in a duct can be related to the displacement of the particle by 


y* + 

Dt at* ax* 


(Bl) 


where an asterisk denotes a dimensional quantity. Assuming a harmonic displacement 
and velocity of the form 


V* ~ I* ~ giw*t* ^g2) 

gives 

v‘ = io)*|* + U* (B3) 

ax* 

where a prime denotes a dimensional quantity with the time function removed. 

With the absence of convection inside the absorbii^ wall, the transverse velocity 
can be expressed as 


V 


♦ 

w 



(B4) 


Again, for harmonic displacements as in equation (B2), the transverse velocity inside 
the wall becomes 


V* = io)*i* 

w ^w 


(B5) 


or 



(B6) 
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The continuity of particle displacement requires that at the interface between the ab- 
sorber and the fluid 


= f»| 

^ I w 

Combining equations (B3), (B6), and (B7) yields 


(B7) 


XT* 9v' 

ia.* 9x* 

The impedance condition at a wall can be written as 

P* 


C = 


3|c 3{c 

P0%^w 


or for harmonic displacements (eq. (B2)) the wall velocity becomes 


(B8) 


(B9) 


V = P 




(BIO) 


Substituting equation (BIO) into equation (B8) yields, for the transverse velocity at the 
wall in the medimn, 


1 /r» , u* 


pjcj? V ico* ax* 


(BID 


In dimensionless form. 


vl =1522 /p ap) 
\ 2nvdxJ 


(B12) 


Substitutii^ equation (B12) into equation (12) in the main body of the report yields the 
displacement boimdary condition at the wall. 


ap 


-i27T7jP 

2M ap 

, iM^ a^p 

3y 

W 


? ax 

w ?2 ;t7J9^2 


(B13) 


w 
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appendix c 


EXIT CONDITION 

The derivation of the exit boundary condition for a soft wall is now presented. The 
goal of the derivation will be to relate the axial pressure gradient at the exit to the im- 
pedance and pressure at the exit. This requires some manipulation of the continuity and 
momentum equations in order to eliminate the acoustic velocities from the boxmdary 
condition. This is necessary since the exit condition is to be used in conjunction with 
the wave equation, which contains only pressure terms. 

At the exit, the impedance relationship can be written 

= -^22? (Cl) 

PqCo u 


or 


C 


(C2) 


At the exit, the momentmn equation (eq. (11)) must also hold; therefore substituting 
equation (11) into equation (C2) yields 


i — + i at the exit plane (C3) 

9x 2jrri 9x 

The axial gradient of u in equation (C3) can be expressed in terms of the transverse 
velocity v by using the continuity equation. Equation (10) can be rewritten as 

— = -i(23T?j)^P - 2nriM - iX (C4) 

9x 9x 9y 

Substituting equation (C4) into equation (C3) yields 


9x 


-i{2nr })( ^ - M j 

\^e ' p ^ M 9v 

1 - 2jt7j(1 - M^) 


(C5) 
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I I IHII II III! I I I 


I II 


The transverse gradient term can also be expressed in terms of pressure. An ex- 
pression for the transverse velocity gradient is obtained by differentiating equation (12) 
with respect to y. 


= i + i _M. 

3y gy2 27T77 3y\ 9x/ 


(C6) 


Substituting equation (C6) into equation (C5) yields 
-i{2Trr])(— - m\ 


3P 

3x 




P + 


iM 


3^P 


iM'^ 




1 - 


27rr7(l - m 2) 0y2 (2 ;t77)^(1 - M^) 


at the exit plane 
(C7) 


The axial gradient in v in equation (C7) can be replaced by the transverse gradient 
in u. Differentiating the axial momentum equation (eq. (11)) with respect to y and 
the transverse momentum equation (eq. (12)) with respect to x and subtractii^ these 
equations yield 


/3u _ 

. iz') 

+ i^( 

^3u _ 

o 

II 

\3y 

dx) 

M ' 

'.3y 

0x/ 


which has a solution 

^ _ Qg(-i277r?/M)x 

3y 3x 


(C8) 


(C9) 


Since the desired solution should also be valid at a Mach number of zero, the constant of 
integration is taken as zero; thus, 


du _ dv 
dy dx 


(CIO) 


This is the condition of irrotationality that is anticipated from irrotational potential flow 
theory. If the potential functions u = 3$/3x and v = 3$/3y had been assumed, equa- 
tion (CIO) could have been assumed immediately without any manipulations. Substitutii^ 
equation (CIO) into equation (C7) yields 
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1? 

3x 


-i{2iT-n)(— - M) 

Ue / p ^ iM 8^P ^ 3^u 

1 - 2nr]{l - M^) 3y^ * (27rr7)^(l - M^) 3y^ 


(Cll) 


The transverse gradient of u can be found from the exit impedance relationship. 
Differentiating equation (C2) with respect to y yields 


3u _ 2 nr] 3P _ 27t? 7P ^^e 

9y ?e 


(C12) 


In this report, is assumed to be constant and independent of y; therefore, 

3u _ 2117] 3P 

9y ?e ^y 


(C13) 


Differentiating one more time gives 


3^u _ 2nv 3^P 
3y^ ?e 3y^ 

Substituting equation (C14) into equation (Cll) yields 


(C14) 


IP 

3x 


e 



1 - 


iM(l + — 


P + 


3/ 3^P| 
2^77(1 - M^) 3 y^ 


(C15) 


For the case under consideration in this report, = 1. Therefore, equation (CIS) re- 
duces to 


^ -i(2?r7?) p ^ iM 3^P 

3x e ~ 1 + M 2;r7j(l - M) gy2 


(C16) 
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APPENDIX D 


FINITE-DIFFERENCE EQUATIONS AND COEFFICIENTS 

The finite-difference equations for the various cells shown in figure 4 can be ob- 
tained by integrating the wave equation over the cell area. 



ay^ 




Kgp + Kg iE^dx dy = 0 


(Dl) 


where the plus sign (+) in the upper limit of integration means to evaluate the param- 
eters along either the upper or right-hand boundary of the integration cell, shown in 
figure 4 by the dashed lines, while the negative sign (-) applies to either the lower or 
left-hand boundary of the integration cell, depending on whether x or y is considered. 

For convenience, equation (Dl) is broken into four separate integrals 


/TS**-// 


K, 


dx dy + / / dx dy + Kc 


+ K, 


/■/■' 

IT 


dx dy 


Each of these integrals can be integrated 


K, 




dx dy 


dx dy = 0 
dx 


(D2) 


+ Kgp 


+ 


x-axis 



(D3) 


The various pressures and their derivatives in equation (D3) are evaluated at the e(%es 
of the cells, designated by the + and - signs. These terms will be different for each 
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cell because the boundary conditions, as well as the size of Ihe cell, vary 


Central Cell 1 (k = 1) 


In this region. 



II 

+ 

+ 

-iE 

9x 

0X 


0X 




iE 

^=iE 


9y 

_ 9y 

3y 







dx = Ax 
dy = Ay 


(D4) 

(D5) 

(D6) 

(D7) 

(D8) 


Substituting equations (D4) to (D6) into equation (D3), performing the required integra- 
tions, multiplyii^ each term by -(Ay/ Ax), and groupii^ the terms into the form given 
by equation (54) yield the following coefficients for the difference equations: 



(D9) 


b 


1 ■ 


-1 


ci = 2 


1 + K 


o- 


2 K2(Ay)" 


d 


1 " 


-1 


(DIO) 

(Dll) 

(D12) 
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(D13) 


e 




- K. 


(Ay)^ 

2 Ax 


Upper-Wall Cell 2 (k = 2) 

In this region, dp/dx]'^, p]"^, and J_ dx remain the same as in cell 1. In this 
region, 3p/9y|^ takes on the value of the boundary condition given by equation (29); 
therefore. 


9y 




w 3x2 


w 



(D14) 


When equation (D14) is substituted into equation (D3), integrating 0p/9xj^ and 
9^p/9x^|^ yields 



3x 


dx = D 


w 


w 





dx=^ 

0X 

w 


w 


+ 




(D15) 


(D16) 


Also, 

2 


(D17) 


Substituting equation (D4), 9p/9y|_ from equation (D5), equation (D7), and equations 
(D14) to (D17) into equation (D3), performing the required integrations, and multiplying 
each term by -2(Ay/Ax) yield the following coefficients: 


ao = a. + W« - 2W, 

^ ^ ^ Ax ,2 


(D18) 



(D19) 
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(D20) 


Cg = Cj - 2Wj Ay + 4Wg 

Ax^ 


d 


2 


= 0 


(D21) 




(D22) 


In equation (54), j takes on the value of m. 


Lower-Wall Cell 3 (k = 3) 

A similar derivation applies to the lower wall, fii performii^ the derivation, re- 
member that the impedance is a vector quantity and must be taken as negative at the 
lower wall if the outward normal of the velocity is not used at the lower wall. The fol- 
lowing coefficients result; 


^3 “ ^2 

(D23) 

bg = 0 

(D24) 


(D25) 

ds = -2 

(D26) 

®3 ~ ®2 

(D27) 


In equation (54), j takes on the value of 1. 


Central Exit Cell Number 4 (k = 4) 


In this region, 3p/ay|t and f_ dy remain the 
9p/0x|‘*" takes on the value of the exit condition given 


same as in cell 1. 
by equation (36); 


In this region, 
therefore. 


iE 

3x 




+ L.i?E 
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(D28) 



0 2 

When equation (D28) is substituted into equation (D3), integrating the 9 p/3y term 
yields 


Also, 



(D29) 




(D30) 



(D31) 


Substituting equations (D5), (D8), and (D28) to (D31) into equation (D3), performing the 
required integrations, and multiplyir^ each term by -2(Ay/Ax) yield the following coef- 
ficients: 


a^ — 2aj 

(D32) 

2K,L„ 

b = -1 - _i_i 

^ Ax 

(D33) 

Kg Ay2 4 KjL2 2K^Lj^ Ay^ 
Ax Ax Ax 

(D34) 

2K^L„ 
d4 = -l--^^ 
^ Ax 

(D35) 

CD 

II 

O 

(D36) 


In the problem, is 1; however, is left in these equations in the event that a differ- 
ent value of might be used at some future date. Also, in equation (54), i takes on 
the value of n. 
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Upper- Corner Exit Cell 5 (k = 5) 


The corner condition combines both the exit and wall impedances in a sii^le cell. 
In this region, 9p/9xlt is again defined by equation (D28). Performing the first inte- 
gration in equation (D3) yields 


K, 


[ 


iE 

9x 


K 


dy = 


l^i- 


U_ 




2 Ax 


2 


+ K1L2 


ie 

9y 


(D37) 


The term 9p/9y|"^ is defined by equation (30) in the body of the report. Unfortunately, 
the second derivative term 9 p/ 9x | ^ in equation (30) can not be evaluated in this cell 
because only two grid points are available in the x-direction. Three grid’points are 
required to evaluate a second-order derivative. 

This problem can be circiunvented by redefining the wall condition in terms of the 
exit impedance. At the exit, for the irrotationality condition (eq. (73)) and the im- 
pedance condition (eq. (C13)), we can write 


9v _ ^ _ 2 nr} 9P 
9y dy 


(D38) 


Substituting equations (D38) and (28) into equation (12) and solving for 9P/9y 5 n.elds 


9y 


wall exit 


i2ff?7P 


M 




iE 

9x 


(D39) 


wall exit 


In, terms of the transformed pressure, defined by equation (18) in the body of this 
report, equation (D39) becomes 


9y 


wall exit 


= N^p 


+ N. 


' 9x 


(D40) 


wall exit 


where 



1 + My 77 J 


(D41) 
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Assuming 



iE 

dx 


wall exit 



(D42) 


(D43) 


then 


iE 

9y 




In addition, p(^ is the same as in cell 4 (eq. (D30)) and 


(D44) 


r 


dx = 


Ax 


(D45) 


/*dy=^ 

” n 


(D46) 


Substituting equations (D30) and (D37) into equation (D3), performii^ the required inte- 
grations, multipl 3 dng each term by -4(Ay/Ax), and grouping the terms into the form 
given by equation (30) yield the followir^ coefficients: 


a- = a. + 4K,L„N9-^+ 2N„-^ 
5 4 1 2 2^^2 2^^ 


(D47) 


bg = 2b4 


Cg = C4 - 2Nj Ay - 2Ng - 4K^L2N^ ^ - 4 KJL 2 N 2 


Ax 


Ax‘ 


dg = 0 


(D48) 

(D49) 

(D50) 
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(D51) 


eg = 0 

In equation (30), j takes on the value of na and i takes on the value of n. 

Lower- Corner Exit CeU 6 (k = 6) 

A similar derivation applied to the lower wall (see discussion of cell 3) yields 


ag - Ug 

(D52) 

bg = 0 

(D53) 

^6 = % 

(D54) 

^6 = ^5 

(D55) 

eg = 0 

(D56) 


In equation (30), j takes on the value of 1 and i takes on the value of n. 
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APPENDIX E 


POLYNOMIAL- EXPONENTIAL INTEGRATION 


The integral in equation (69) is expressed as the sum of the integrals between each 
grid point: 



p(s, y)e^^® ds 


rx.+Ax 

=y J p(s, y)e^^® ds 


(El) 


A third-order pol 3 momial fit of the pressure p is constructed from the solution of the 
wave equation 


2 3 

p(s, y) = CKq + Q!l(s - Xq) + 0-2(3 - Xq) + o-3(s - Xq) 


where 

Starting point: 


Central points: 


Xq = 0 if X. = 0 


Xq = Xj - Ax if 0 < X. < Xg - Ax 


(E2) 


(E3) 


(E4) 


End point: 


Xq = Xj - 2 Ax if X. = Xg - Ax 


(E5) 


Substituting equation (E2) into equation (El) gives 

n 


/ Xq rxj+Ax 

p(s, y)e^^® ds = \ / [oq + q-j(s - Xq) + 0-2(8 - Xq)^ + 0-3(3 - XQ)^]e^^® ds 

4^ (E6) 


The integral of the polynomial in equation (E6) will be integrated exactly by usii^ 
formulas given in standard integration tables (ref. 16). However, to simplify the inte- 
gration, let 
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Z = S ~ Xf 


(E7) 


Therefore, equation (E6) becomes 

n 


/ 


"Xi +Ax-Xfi 


^ _ . 3 ^ 

i 


|ofo + dz (E8) 

Furthermore, to ease the numerical integration, let 


X. - Xq = N Ax 


X. + Ax - Xq = (N + l)Ax 


where 

Starting point: 


Central points: 


End point: 


N = 0 if X. = 0 


N=1 ifO<x. <x - Ax 
1 e 


N = 2 if X. = x^ - Ax 
1 e 


Therefore, equation (E8) becomes 

n 




p(s, y)e^^® ds 


(E9) 

(ElO) 


(Ell) 


(E12) 


(E13) 


(“0 “1=^ * “2“^ * o, 32 ’)e‘’’=’ dz 


(E14) 


By usii^ equations (565. 1, 567. 1, 567. 2, and 567. 3) from reference 16, the integral 
for equation (E14) can be expressed as 
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n 

® p(s, ds =^ e - l) + o,^ Ax[e^‘i^ ^^(N + 1 - e) - (n - c)] 

i 

+ «2 ^^[(N + 1)^ - 2e(N + 1) + 2e^] - (N^ - 2eN + 2€^)} 

+ 0-3 Ax^p*;^ ^^[(N+ 1)^ - 3e(N+ 1)^+ 6€^(N+ 1) - 6e^] 


(N^ - 3eN^ + 6e^N - Ge^)] 


(E15) 


where 


e = 


i(p Ax 


(E16) 


The expressions for the coefficients o-q, a>p o-g, and otg can be evaluated directly 
from the formulas 


“0 = % 

Q! 2 = 3 Qj - 3Q2 + Q3 
8Q2 “ “ ^^3 


Q!o = 


2 Ax 


"3 = 


Qg " 2Q2 + Qj 


(E17) 

(E18) 

(E19) 

(E20) 


2 Ax^ 


where 


Ql = 


Qq = p(xq» y) 

p(Xq + Ax, y) - p(xq, y) 


Ax 


(E21) 

(E22) 
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p(xq + 2 Ax, y) - p(xq, y) 

(E23) 

2 Ax 


p(xq + 3 Ax, y) - p(xq, y) 

(E24) 

3 Ax 
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appendix f 


SAMPLE PROBLEMS 


The appendix contains complete pressure and velocity printouts for four sample 
problems. These printouts can be used to debug and check computer programs. 


Problem 1 


Plane 

pressure 

wave 



(a) 


The input is 


T] = 0.6 

tf' = 0 (wave envelope concept not used) 


M= 0 


L*/H* = 0. 5 

= 0. 16 - i 0. 34 
P(0, y) = 1 


n = 5 


m = 10 
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The calculated values are 


Ax = 0. 1 
Ay = 0. 1111 
Eq = 1. 30 

= 0. 128 (x = 0. 5) 

A 

AdB = -10. 04432 


k 




bk 





<*k 



1 

-.123*01 

.cno 

-.100+01 


.000 

.429+01 

.noo 

-.100+01 


.000 

-.123+01 

2 

-. 123+01 

.030 

-.200+01 


.000 

.224+01 

.949+00 

.000 


.000 

-.123+01 

5 

-.123+01 

.000 

.CQC 


.000 

.22R+D1 

,949+00 

-.200+01 


.uoo 

-.123+01 

4 

-.247+01 

.000 

-.100+01 


.000 

.429+Cl 

.931+00 

-.100+01 


.COG 

.000 

5 

-.247+01 

.000 

-.200+01 


.000 

.228+01 

.188+01 

.UOG 


.bOO 

.OOP 

6 

-.24 7+01 

.000 

.000 


.000 

.228+01 

.188+01 

-.200+01 


.000 

.noo 


i 

X position 

3 

|P| 

j |Pl 

j 

|P| 

j 

|Pl 

j 

IPI 

1 

.0000 

1 

.100+01 

2 ,100+01 

3 

.100+01 

4 

.100+01 

5 

.lPC+01 

2 

.lUOO+OD 

1 

.108+01 

2 .725+00 

3 

,793+00 

4 

.880+00 

5 

.926+00 

3 

.ZuDO+OO 

1 

.m5+ci 

2 ,518+00 

3 

.578+00 

4 

.717+00 

5 

.792+DQ 

4 

.3000+00 

1 

.793+00 

2 .318+00 

3 

.393+00 

4 

.S6C+0P 

5 

.650+00 

5 

.4000+00 

1 

.599+00 

2 .169+00 

3 

.262+00 

4 

.446+or 

S 

.541+00 

6 

.5000+00 

1 

.598+00 

2 .144+00 

3 

.213+00 

4 

.404+00 

5 

.500+00 

1 

.0000 

6 

.100+01 

7 ,100+01 

6 

,no+oi 

9 

. 100+01 

10 

.lDQ+01 

2 

. lGOO+00 

6 

.926+OC 

7 ,88C+0D 

8 

.793+00 

9 

.725+00 

10 

.108+01 

3 

.2000+00 

6 

.792+00 

7 .717+00 

8 

.578+00 

9 

.518+00 

10 

.105+01 

4 

. 3CC0+OO 

6 

♦650+00 

7 .560+00 

8 

.393+00 

9 

.318+or 

10 

. 703 +OL 

5 

.4000+00 

6 

.541+00 

7 .446+00 

8 

.262+00 

9 

.169+00 

10 

.599+00 

6 

.5000+00 

6 

.500+00 

7 .4^4+00 

8 

.213+nn 

o 

.144+00 

10 

.598+00 


i 

X position 

j 

p(l) 

3 

p(l) 

3 

p{l) 

j 

p(l) 

j 

p(l) 

1 

.0000 

1 

.1<"C+01 

7 

.l^O+Ol 

3 

,100+01 

4 

,100+Cl 

5 

.100+01 

t. 

.ILPO+QO 

1 

.341+00 

7 

,574+00 

3 

,752+0" 

4 

,860+C" 

5 

. 9P9+G0 

3 

.2CjCQ+C0 

1 

-.511+OC 

7 

.111+00 

3 

.454+0"' 

4 

,6fafa+un 

5 

.729+00 

4 

. 3000+00 

1 

-.758+50 

? 

-.142+00 

3 

,216+OC 

4 

.421+00 

5 

.515+00 

5 

.4000+00 

1 

477+DU 

2 

-.166+30 

3 

,725-01 

fa 

.228+C" 

5 

.303+00 

6 

.5000+00 

1 

-,131+OC 

? 

-.109+03 

3 

-,142-01 

fa 

.670-01 

5 

,110+00 

1 

.COOO 

6 

.iro+ci 

7 

.If'O+Ql 

8 

,irc+oi 

9 

.100+01 

10 

.IPC+Ol 

2 

. IbOO + UO 

6 

.9"9 + 0L 

7 

,660+00 

8 

,752+On 

Q 

.574 + Cf' 

IP 

.341+Qb 

3 

.2000+00 

6 

.729+DO 

7 

.644+00 

8 

.454+00 

9 

.111+0'^ 

10 

-.511+00 

4 

.3000+00 

6 

.515 + CC 

7 

,421+00 

S 

,216+00 

9 

-.142+or 

10 

-.758+00 

5 

.4LD0+00 

6 

.3P3+CP 

7 

.226+00 

8 

,725-01 

9 

-.166+00 

10 

-.477+00 

6 

.sloo+co . 

6 

.110+00 

7 

.670-01 

8 

-,142-Pl 

0 

-.109+cn 

10 

-.131+00 


i 

X position 

j 

p(2) 

j 

p(2) 

j 

p(2) 

j 

p(2) 

j 

p(2) 

1 

.0000 

1 

.OPC 

2 

.ono 

3 

.ono 

4 

.COD 

5 

.QPD 

2 

.1003+00 

1 

-.103*01 

2 

-. 443 +DO 

3 

-.250+rn 

4 

-,190+or 

5 

-.174+00 

3 

.200CI+UO 

1 

-.916*00 

2 

-.5^6+00 

3 

-,358+rn 

4 

-,316+cn 

5 

-,31C*C0 

4 

.3000*00 

1 

-.234+CO 

2 

-.285*00 

3 

-,327+CO 

4 

-.37D+CP 

5 

-.396+cG 

5 

,4bC0+D0 

1 

,363+00 

2 

-.289-01 

3 

-.251+On 

4 

-.384+0" 

5 

-.443+Ou 

6 

.50C0+G0 

1 

. 544+03 

2 

.935-31 

3 

-,212 + P'^ 

4 

-. 398+cr 

5 

-.4P8+0P 

i 

,CCOO 

6 

.bPP 

T 

.oro 

a 

• ^rc 

0 


10 

.bfu 

2 

,1000+00 

6 

-. 174+n'-- 

7 

-.140*30 

8 

-.250+02 

0 

- .44 3 + cr 

10 

-.133+01 

3 

.2000+00 

6 

-. 31C + 0t- 

7 

-,316+on 

6 

-,356+P" 

9 

-,5"6+jr 

IP 

-.916+30 

4 

.30P0+C0 

6 

-. 306 +C 3 

7 

-.37U+C2 

6 

-.327+0" 

0 

-.285 + 'l" 

13 

-.234+00 

5 

. 4C0O+OO 

£ 

-,448 + C-D 

7 

- ,384+Cn 

8 

-. 29 1 ♦ ri" 

9 

-.249-31 

10 

.365*00 

6 

.9003*00 

6 

-.4P8 + C L 


-.346+32 

e 

-,212 + f'" 

° 

.935-01 

1C 

.5 P4 4!;p 


.oon 
.000 
.000 
• 000 
• oon 
.000 
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i 

X position 

3 


3 

u(l) 

j 

u<« 

3 

,u<« 

3 

. u(« 

1 

*'cno 

1 

.177+L? 

? 

.b?7+01 

3 

.302+01 

4 

.2f’7 + Ql 

5 

.i« 2 +ni 

2 

, ICCJ+uO 

1 

B?+a 1 

2 

.3^9+01 

3 

.209+01 

4 

.172 + 01 

5 

.163+01 

J 

. 2lOU+uO 

1 

- . 4 + u 1 

2 

-.i^a+ji 

3 

.3'lb + OO 

4 

. BBl +00 

5 

.111+01 

k 

. 2Ln0*00 

1 

-,7ie*L 1 

2 

-.4.73+01 

3 

-.b77+m 

4 

.275+br 

5 

.6'^9+On 

S 


1 

-.31 5> + ai 

? 

-.1'1+Ql 

3 

-.472+00 

4 

. IfrO+OD 

S 

.455+00 

b 

.£lDJ+aO 

1 

“.2i‘>6*n0 

2 

-.2^6+00 

3 

.545-01 

4 

.301+00 

5 

.440+00 

1 

.CLOD 

6 

. 162 + Cl 

7 

.207+01 

a 

.3'’? + ni 

9 

.627+01 

10 

.177+02 

2 

. lCOO+00 

6 

.lb3*01 

7 

.172+01 

8 

.209+01 

o 

.309+01 

10 

.482+01 

i 

.2LO0+00 

6 

.11 1*01 

7 

,a«l+00 

B 

. 3Ofe.P0 

9 

-.108+01 

10 

-.499+01 

>t 

. 3LDO+CO 

6 

• 6 “1 9 +00 

7 

.275+00 

8 

-.b77+no 

9 

-.273+01 

10 

-.718+01 

5 

.‘*lOJ*QD 

6 

.45S*L0 

7 

. lbO + 30 

8 

-.472+00 

9 

-.151+01 

10 

-.315+01 

6 

. 5CaC+00 

& 

.440+00 

7 

.301+00 

8 

.545-01 

o 

-.2>^6 + lj0 

10 

-.266+00 


1 

X position 

3 

CM 

3 

u(2) 

3 

U® 

3 

u(2) 

3 

u(2) 

1 

.0000 

1 

-.795+CO 

2 

-.268*01 

3 

-.168+01 


-.724+00 

5 

-.239*00 

2 

.1000+00 

1 

-.108*02 

2 

-.570*01 

3 

-.299+01 

4 

-.191+01 

5 

-.145+01 

3 

.2000+00 

1 

-.609+01 

2 

-.380+01 

3 

-.280+01 

4 

-.228+01 

5 

-.205+01 

4 

. 3LD0+00 

1 

.620+00 

2 

-.128+01 

3 

-.191+01 

4 

-.211+01 

5 

-.217*01 

5 

.4000+00 

1 

.302 + 01 

2 

.393*00 

3 

-.105+01 

4 

-.174+01 

5 

-.202+01 

6 

.500Q+00 

1 

.127+01 

2 

.744-01 

3 

-•679+00 

4 

-.152*01 

5 

-.185+01 

1 

.0000 

6 

-.239+00 

7 

-.724+00 

6 

-.168+01 

9 

-.268+01 

10 

-.795+00 

2 

.1000+00 

6 

-.145+01 

7 

-.191+01 

8 

-.299+01 

9 

-.520+01 

10 

-.108+02 

3 

.2U00+DG 

6 

-.2f'5+01 

7 

-.228+01 

8 

-.280+01 

9 

-.380+01 

10 

-.609+01 

4 

. 3CD0+00 

6 

-.217+01 

7 

-.211+01 

8 

-.191+01 

9 

-.128+01 

10 

.620+00 

5 

.4000*00 

6 

-.202+01 

7 

-.174+01 

3 

-.105+01 

9 

.393+00 

10 

.302+01 

6 

.5000*00 

6 

-.185+01 

7 

-.152*01 

8 

-,879+CO 

9 

.744-01 

10 

.127+01 


1 

X position 

3 

yd) 

3 

yU) 

3 

yd) 

3 

yd) 

3 

yd) 

1 

.0000 

1 

.ono 

2 

.oou 

3 

. jno 

4 

.000 

5 

.000 

2 

.1000+00 

1 

-.820+01 

2 

-.273+01 

3 

-.681+00 

4 

-.258+00 

5 

-.507-01 

3 

.2000+00 

1 

-.552+01 

? 

-.213+01 

3 

-.684+00 

4 

-.141+00 

5 

-.898-02 

4 

.3000+00 

1 

,560+00 

2 

, 332+00 

3 

.3P9+00 

4 

.328+00 

5 

.126+00 

S 

.4000+00 

1 

.466+01 

2 

,251+01 

3 

.152 + 01 

4 

,864+aP 

5 

.285+00 

6 

.SLOO+DO 

1 

.553+01 

2 

.340+01 

3 

.215 + 01 

4 

.122+01 

5 

.397+00 

1 

.0000 

6 

.000 

7 

.000 

8 

.uOO 

9 

.000 

10 

.000 

2 

. lDOD+00 

6 

.507-01 

7 

.258+00 

8 

.681+00 

9 

.273+01 

10 

.820+01 

3 

.2C0D+DD 

6 

,898-02 

7 

,141+00 

8 

,684 + nf^ 

9 

.213+01 

10 

.552+01 

4 

. 3000+00 

6 

-.1 26+00 

7 

-, 328+30 

B 

-.389+00 

9 

-.332+30 

ID 

-.560+00 

5 

.4000+00 

6 

-.285+00 

7 

-.664 +00 

8 

-.1 52+01 

9 

-.251+01 

10 

-.466+01 

6 

.6L0D+00 

6 

-.397+00 

7 

-,1?2+01 

8 

-.215+01 

9 

-.340+01 

10 

-,553+Di 


i 

X position 

3 

y(2) 

3 

y(2) 

3 

y(2) 

3 

y(2) 

3 

y(2) 

1 

.0000 

1 

.000 

2 

.000 

3 

.JOG 

4 

.000 

5 

.000 

2 

.1000+00 

1 

,224+01 

2 

.197+01 

3 

,129+01 

4 

.691+00 

5 

.217+DQ 

3 

.2000+00 

1 

.739 + 01 

2 

.400+01 

3 

.227+01 

4 

.119+01 

5 

.370+00 

4 

.3000+00 

1 

.715+01 

2 

.411+01 

3 

.242 + 01 

4 

.130+01 

5 

.410+00 

5 

. 4000+00 

1 

.304+01 

2 

.253+01 

3 

.178 + 01 

4 

.103+01 

5 

.337+00 

6 

.5000+00 

1 

-.538+00 

2 

,799+00 

3 

,876+00 

4 

,581+GC 

5 

.197+00 

1 

.0000 

e 

.OPQ 

7 

.uOO 

8 

.000 

9 

,000 

10 

.000 

2 

.1000+00 

6 

-.217+00 

7 

-.691*00 

8 

-.129+01 

9 

-.197+01 

10 

-.224+01 

3 

.2000+00 

6 

-.370+00 

7 

-.119+01 

8 

-.227+01 

9 

-.400+01 

ID 

-.739+01 

4 

. 3000+00 

6 

-.410+00 

7 

-.130+01 

8 

-.242+01 

9 

-.411+01 

10 

-.715*01 

5 

.4000+00 

6 

-.337+00 

7 

-.103+01 

8 

-.178+01 

9 

-.253+01 

10 

-.304+01 

6 

.5000*00 

6 

-,197+OD 

7 

-.5»l+00 

8 

-.876+00 

9 

-.799+00 

10 

.538+00 


Problem 2 
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Problem 2 is the same as problem 1 except that the wave envelope concept is used. 
The input is 


77 = 0. 6 
if = 0.6 
M = 0 

L*/H* = 0. 5 

= 0. 16 - i 0. 34 

p(0, y) = 1 

n = 5 
m = 10 


The calculated values are 


Ax = 0. 1 


Ay = 0, 1111 


Eq = 1.21 

= 0. 153 (x = 0. 5) 
AdB = -9. 01164 


k 




*>k 








1 

-.123*01 

- .465*00 

-.100*01 


.000 

.447*01 

.non 

-. ino*oi 

.LQC 

-.123*01 

.465*00 

2 

123*01 

-.465*03 

-.200*01 


.oco 

.745*01 

, Q 14 9 ♦ n n 

.uOC 

.LOC 

-.123+01 

.465*00 

5 

-.123*01 

-.465*00 

.COO 


.DOD 

.745+01 

,g49*0D 

-.200*01 

.0"C 

-.123*01 

.465*00 

4 

-.247*01 

-.931*00 

-.100*01 


.000 

.447*01 

.031*^0 

-.100*01 

.000 

.roc 

.000 

5 

-.247*01 

-.931*00 

-. 2 ro*Ll 


.000 

. 745*01 

. 184 + 01 

.COO 

• bOC 

.COC 

.000 

b 

-.247*01 

-.931*00 

.□□0 


.000 

.745*01 

.184*01 

-.200*01 

.one 

.nor 

.TOO 


46 



i 

X position 

j 

|P| 

j 

IP| 

j 

I^I 

j 

IP| 

j 

IP| 

1 

.LUOO 

1 

.ino+oi 

2 

.l''D + 01 

3 

. 1"C+01 

u 

. 10(? + J1 

5 

. 1 '’C + 1- I 


. Uj03+00 

1 

.1 in + Di 

2 

.738+on 

3 

♦ 799+P" 

u 

. 3 4 0 + (J n 

5 

. 9?3 + uL 

l 

. ?COD*CiO 

1 

.1C5+01 

7 

.5-'3*0Q 

» 

, SoC+"3 

u 

. 722*.." 

5 

. 793+ ’l 

<4 

, 30rG*u0 

1 

.7 98+or: 

J> 

.33C+OC 

■» 

.4'*6 + C n 

u 

,S7r ♦ ;P 

5 

. b 7 + >. 

S 

. 4bDD-»ljO 

1 

.635+CC 

? 

. 167+3" 

3 

,273+n" 

4 

.449+ 

5 


6 

.5L00*Cn 

1 

.6«1+CG 

2 

.167+30 

3 

.224 + n 

4 

. 4 1 8 + '. r 

5 

.sib+"r 

1 

.^LOO 

e 

.I'^n + oi 

7 

.l^C+01 

3 

.i('0 + 01 

9 

. 1"0+^1 

I'J 

. I ''0 + ., 1 

Z 

. icoo*un 

6 

.933+OC 

7 

.aPC+00 

5 

.799+0" 

o 

.736+ 

n 

. 1 1 ♦ M 

3 

.2CG0*Q0 

6> 

.793+LC 

7 

.7'*2 + 3T 

6 

, b'SC ♦ 

« 

.473+ r 

13 

. 1 ' 5 + ' i 

«4 

. luOO +DO 

6 

, 6 7 ♦ C Ci 

7 

.570+30 

8 

,H''6 + p" 

n 

. 3 3r+ ;r 

n 

. 7''P+ 'L 

5 

.‘*cCIO*LO 

6 

.5«=M*ur 

7 

.4'-9 + jn 

4 

.273*''- 

Q 

.167+." 

10 

,i7^* . 

& 

. SL3C ♦•iJC 

b 

. 5 1 6 + G C 

7 

.418+30 

a 

.274+''» 

9 

. 167 + C " 

IC 

,c'^l+ ■ 


i 

X position 

i 

p(« 

j 

p'« 

3 

p(« 

j 

p(« 

i 

p<« 

1 

. 30DD 

1 

.100+01 

2 

,iro+Di 

3 

. Kc+n 

4 

. no+bi 

5 

. 1 "u + 01 

2 

. U DO+OD 

1 

.779+00 

7 

.712 + 30 

3 

.798+0" 

4 

.370 + .?" 

5 

. V ''6 + 4C 

3 

. 2L0Q + U0 

1 

.334+GO 

7 

.449+30 

3 

,5®8+rr 

4 

.694 + c.r 

5 

.7'-l+ ' 

4 

. 3CC0 + C0 

1 

“.860-01 

2 

,212+30 

3 

.403+Cr 

U 

,529+rr 

5 

, 5944 ' " 

5 

.4000+00 

1 

“. 4"4 + 00 

2 

,236-G1 

3 

.270+0" 

4 

.4 1 7+C P 

4 

.4 48 + .).. 

& 

. EOOO + uO 

1 

-.541+00 

? 

-.531-01 

3 

,221+on 

4 

. 775 + .:r 

5 

.4«-3 + "f' 

1 

.COCO 

6 

. 1"0 + 01 

7 

,100 + 01 

8 

. 1"0 + P1 

o 

. inn + .ji 

10 

, i"C: + r 1 

2 

. ILOO + OO 

t 

.9"8 + 0a 

7 

.d‘"L+DO 


. 798 ♦ 

Q 

.712+J" 

in 

. ?r9 + '-L 

3 

. ?uOo + on 

6 

.751+00 

7 

.694+30 

8 

,S“8+r" 

9 

.449 + . 1’ 

10 

,3''4 + . 0 

4 

.3003+00 

6 

. 5 9 4 + 0 a 

7 

.579+30 

6 

.4"3+00 

9 

.21 2 + ..r 

10 

-.560-31 

5 

.40CC+O0 

6 

.488+00 

7 

.417+00 

3 

.270+0" 

9 

.236-01 

10 

-.4"4 + .)L 

b 

.5000+00 

6 

.453+00 

7 

.378 + 00 

8 

.221+0" 

e 

-.•^’l- .1 

1 c 

-. 54 1 + M. 




J 

d^2) 


0^2) 

j 

d<2) 


„(2) 


J2) 

i 

X position 

P 

j 

P 

P 

j 

P 

i 

P 

1 

. i.D"0 

1 

.uGO 

2 

.u"C 

3 

..-"C 

4 


5 

. j" . 

2 

, iL"a + u3 

1 

-. e ip+co 

? 

-.196*30 

3 

,413-"1 

4 

.1^2 + 1 ' 

5 

. m+ ;> 

3 

. 2 l.33 +U0 

1 

-. in + 01 

? 

-. .147 + 30 

3 

,44^-01 

4 

.19“+." 

s 

.?"'7 + . 

4 

, :l"G+ji: 

1 

-.793+00 

7 

-.293*00 

T 

,4»4-r 1 

4 

.PlC+wT 

6 

.2^ '.+ ' 

5 

,4Lcn + (j0 

1 

- . 4 9 P ♦ C 0 

7 

-.146+30 

3 

.44?-"! 

4 

.19 1+ " 

5 

. 2 " 2 ♦ 

6 

. 5i.nD + or 

1 

36 1 +0C 

2 

149+JP 

3 

. 39P-L 1 

4 

.173,+ ' 

5 

.7'*7 + - . 

1 

. lU.3 0 

6 

,u"C 

7 

.L "O 

4 



,■'1 r 

1 C. 


2 

. iLor +UC 

6 

.167+r.C 

7 

,1 ’2+3': 

8 

.4 1 7-n 

T 

-.196+ ■" 

ir 

-.PI P+"i 


.runo+uc 

6 

.^n+LO 

7 

. 19fa+00 

6 

.4‘=9-Ll 

Q 

- . ?67 + , " 

1 - 

-. 1 ' ! + ;• 1 

4 

. 5c.oo + on 

6 

... “C + E.3 

7 

.2’L+C3 

8 

.4 94-" 1 

0 

-.743+ r. 

i ' 

-.700+ • 

5 

, 4 L 00 + 60 

6 

.^62*00 

7 

.191+0'' 

P 

. 4 4 3-r 1 

9 

- . 1 “6+ . 

1 n 

-.49. + 

6 

. + Lin 

6 

.297 + 1 c 

7 

. I’E +0" 

6 

. i9"-u 1 


-.159+ , " 

i 0 

-. 3M +Pi 


1 

X position 

j 


i 


i 


j 


j 


1 

. 0L"O 

1 

. 1 96 + 62 

7 

.599 + 01 

5 

..;96 + r 1 

u 

.2"9 + ‘, 1 

5 

. m + - 1 

2 

. I 1 .OO + JO 

1 

.749+61 

3 

. 499+01 

5 

.299+01 

4 

.232+ Jl 

5 

. 7'^ ^ + f J 

3 

.2oro+uo 

1 

. 5 7 4 ♦ u 0 

? 

.192+31 

3 

.2 71 + f‘l 

4 

. 7?6+. 1 

5 

. 7 7 ♦ ,1 

4 

. 30nG+LO 

1 

-.340+01 

2 

, 1 9 3 + 3 P 

3 

.153+01 

4 

.7''2 + l 1 

5 

. 4 7 'i ♦ 1 

5 

. 4L00 + 0O 

1 

-.319+01 

? 

-.191+03 

3 

.171+31 

4 

.181+^1 


,2' 6+31 

6 

.5i_nc + t.n 

1 

-.172+01 

2 

-.76u-Dl 

3 

. f 4 9 ♦ 0 " 

4 

.138+61 

0 

.164 + 01 

1 

.1)1,03 

6 

. 1P7 + C1 

7 

.2r9+Dl 

8 

.296+PI 

0 

. 5 8 9 + :. 1 

13 

.19 6 + ' : 


. 1C Dc + on 

6 

.2"8+Gl 

7 

.2’2 + 01 

6 

.299+Pl 

9 

.4“ 9+ ,1 

1 6 

.769 + r'l 

3 

. 2uDQ+60 

6 

.227+01 

7 

...76+ni 

5 

. 4 2 1 ♦ r 1 

9 

. 19 2 + 1. 1 

i L 

.S74 + -1 

4 

. 3600+60 

6 

.278+01 

7 

.^"6+01 

b 

.1'^ 3+Pl 

9 

.15 3+ r 

lU 

346+'-'l 

5 

.4Lru+cn 

6 

.2"6+01 

7 

.191+61 

8 

.171+01 

0 

-.191+." 

1 6 

- . 369+P 1 

6 

.SLOO+uO 

6 

.164 + 61 

7 

.I’E+Dl 

8 

. 649 + r ■’ 


- .76' -,'l 

IJ 

- . 1 72 + ' 1 


i 

X position 

j 


j 


j 


i 


j 


1 

. jlP J 

1 

- . 9 7 4 + L C’ 

7 

-. i‘'5 + Jl 

3 

- . 1 *’ " + M 

4 

-.‘-3 3+ ' 

6 

-.316+ 1 

2 

. 1630+40 

1 

-.756+0 1 

7 

-.3' 2+Jl 

3 

-.i97+r 1 

4 

-.113 + 1 1 

5 

- . 7 ’ 2 ♦ 1 

3 

. 26"3+L3 

1 

-. 8 '7 + L 1 

7 

-.jbu+Cl 

3 

- . 1 9 9 + r 1 

4 

-. 1 : 4+ .1 

5 

- , 7 1 

4 

, 3oDO+aO 

1 

-. 666+01 

2 

-.315+61 

3 

14 6 + 1-1 

4 

-.64 7+ " 

5 

-.312+ . 

5 

. 4636+60 

1 

-.40?+61 

? 

-.271+01 

3 

-.ttl 1 + ''" 

4 

-.414-.,] 

5 

. 3" V. + '.) 

£, 

. 5600 +00 

1 

-, 178+r 1 

2 

-.4 06+ 17, 

3 

. 1 9 f • ♦ f '' 

4 


5 

.V 1 P ♦ . 

1 

.063: 

6 

-.3 16+6 a 

7 

-.fe'*5 + 6 P 

8 

- . 1 p a ♦ r 1 

9 

- . 3 " 5 + : 1 

1 6 

- , 9 " 4 + ’ i. 

2 

. KOC + uO 

6 

7 ’2 + ( 6 

7 

-.113+01 

F 

197+Pl 

9 

- . 3 = 2 + Gl 

1 G 

-.7=6+ M 

3 

. 2600 +uD 

6 

- . 7 " 0 ♦ u D 

7 

-,1"4+01 

P 

1 85 +r 1 

9 

- . ’b" ♦: 1 

10 

-.P"7 + f 1 

4 

. 30no+60 

6 

3 12+OP 

7 

- .b47+cn 

S 

-.146 + n 

0 

-.315+61 

IP 

-.666+31 

S 

. 4 600+60 

6 

. irp + cr 

7 

-.415-01 

8 

-.311+"" 

9 

-.27 1 + -jl 

13 

-. 49 ; .pj 

6 

.5600+63 

6 

,9 18+60 

7 

,fa7u+:0 

8 

.1=6+3" 

u 

- .406+ 'P 

IP 

-. 1 7e + r.i 


i 

X position 

3 

v<« 

i 

v<« 

j 

vl« 

3 


3 

v'« 

1 

.0630 

1 

,6PD 

2 

.630 

3 

.CPC 

4 

.POO 

5 

.JPG 

2 

. 1600+60 

1 

-.639+01 

2 

-.316+01 

3 

-.123+01 

4 

- .46 1+GP 

5 

1 17+cr 

s 

.2600+00 

1 

-.912+01 

7 

-.420+01 

3 

-.196+01 

4 

- . P 5 8 * 0 r 

5 

-.241 +G0 

4 

. 3uDU + 00 

1 

-.634+61 

? 

-. 395+01 

3 

-.197+01 

4 

- ,092+6" 

5 

-.331+60 

5 

.46Da*on 

1 

3"0+6l 

? 

“,247+01 

3 

-.17C+C1 

4 

- . 9 7 P + 0 P 

5 

-.513+00 

6 

. 5600+00 

1 

-.158+61 

2 

199 + 01 

3 

-.196+01 

u 

-.945+uP 

5 

-.312+00 

1 

.COCO 

6 

,600 

7 

.DPD 

8 

.COO 

9 

.000 

10 

, "■"0 

2 

.1600+00 

6 

. 117 + 00 

7 

.461+00 

6 

.123+01 

0 

,318+61 

1C 

. 339 + 01 

3 

.2uDQ+aO 

6 

,241+00 

7 

,898 + On 

a 

.196+01 

9 

.420+01 

10 

.912+01 

4 

. 3LOO + OD 

6 

.3P1 + 00 

7 

,992+00 

a 

.197+n 

9 

.355+61 

10 

.634+01 

5 

,4600+00 

6 

.313 + 00 

7 

,970*00 

8 

.17C + P1 

9 

. 247+01 

IP 

. 330 + 01 

6 

.5600+00 

6 

.312 + 00 

7 

,945*00 

8 

. 196 + 01 

9 

. 199+01 

10 

. 153 + 31 



1 

X position 

j 

<1 

j 

y(2) 

J 

v(2) 

3 

''a' 

j 


i 

.cono 

1 

.one 

? 


3 

.wHC 

4 


5 

.IJOO 

2 

. icon^LH 

1 

1 la+ni 

? 

.719+3'’ 

3 

.816 + f’*' 

U 

.5 17 + C'i 

5 

.173+00 

3 

.2O0D+U0 

1 

.122+Jl 

? 

.173+01 

3 

.114+01 

4 

.743+cr 

S 

.258+Oc 

<4 

.3L0n+UQ 

1 

.379+ul 

2 

.2''2+31 

3 

.175+01 

4 

.862 + ^'’ 

5 

,204 + '\ 

5 

.4000+00 

1 

.505+01 

2 

.277+31 

3 

. 185+01 

4 

.962+nr' 

5 

.317+00 

6 

,i>G00+CD 

1 

.50" + Cl 

2 

.314*31 

3 

.1S4+P1 

U 

.ini+Gl 

5 

. 328+fi 

1 

. ocoo 

6 

.u’^D 

7 

.Of'C 

6 


n 

.one 

10 

,000 

2 

. ILOO + QO 

fe 

-.173+00 

7 

-.bl7+or. 

a 

-.G l6 + 0"j 

o 

-.719 + 01' 

10 

.118+01 

3 

.2u00+CD 

6 

-.2'-a + oc 

7 

-.746+OC 

a 

-.114+01 

9 

-.133+01 

10 

-.172+01 

4 

,3lD0+U0 

& 

-.294+OC 

7 

-,b62+Ln 

a 

-.139+01 

9 

-.202+01 

IC 

-.379+01 

5 

.4DD0+0C 

6 

-.317*00 

7 

-.9ft2+Ejr 

8 

-.189+-1 

o 

-.277+lJl 

10 

-.5'5+Dl 

6 

.suno+uo 

6 

-.5?B+f u 

7 

-.1^1+01 

3 

- , 1 « 4 + : 1 

o 

-. 314+ul 

10 

-.5«=3 + 01 


Problem 3 


Plane 

pressure 

wave 


.’T'' 

1 i 


Wi'u//// 


---t- 


1=0 {use ?„ = «} 


(c) 


The input is 


?7 = 0 . 6 

rf' = 0 (wave envelope concept not used) 
M = 0. 5 

= 0. 071 - 0. 151 (upper) 

= 1x10^® (centerline) 

P(0, y) = 1 
n = 10 


m = 10 
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The calculated values are 


Ax = 0. 05 
Ay = 0= 05556 
Eq = 2. 64 

Ejj = 0. 725 (x = 0. 5) 
AdB = -5. 60602 


k 










®k 

1 

926*00 

- .116*00 

- .100+01 

.000 

. ■'ai + Di 

, ODD 

-.1^0*01 

.Ono 

-.926 +rC 

. 116*00 

2 

. 122*02 

-.137*02 

-.200*01 

.000 

-.704*02 

,161*^2 

. LlOu 

. lOC 

. 179*02 

-.137*01 

3 

-.926*00 

-.116*00 

.000 

.000 

.781+01 

.316-29 

-.200*01 

. L'^O 

-.926+00 

. U 6 + 00 

4 

185*01 

-.233*00 

-. ICD+ul 

-.796* u 1 

.781*01 

.164*02 

-.100*01 

-.796*01 

.<^□0 

.000 

5 

.282*02 

-.193*02 

-.200*01 

-.159+ u ? 

-.335*02 

.741*02 

. JOG 

.000 

.000 

.000 

6 

-.165*01 

-.233*00 

.000 

.000 

. 381*01 

.164*02 

-.200*01 

- .159*02 

."<00 

.000 


i 

X position 

j 

|5| 

i 

IP| 

i 

1^1 

j 

|Pl 

3 

m 

1 

.GGOO 

1 

. 1 no*oi 

7 

.1^0+01 

3 

, 100*01 

4 

, lOO+Ol 

5 

. l'’0 + 01 

2 

.5000-ul 

1 

. 9‘^2+GO 

2 

.951*00 

3 

.949*0'^ 

4 

.945+DC 

5 

,941*CC 

3 

. lGDD + 00 

1 

.657+OC 

2 

,856*00 

3 

,S52*nn 

4 

. 847+00 

5 

.843*00 

4 

.1500*00 

1 

.725*00 

7 

.723*00 

3 

.719*00 

4 

.715*00 

5 

.715*00 

5 

.2000*00 

1 

, 567*00 

7 

.566*00 

5 

.563*00 

H 

.563*0,0 

5 

.571*00 

6 

.2500*00 

1 

.urs+oo 

2 

,4^5+On 

3 

.406+0^ 

4 

,414* or 

5 

.433+00 

7 

. 3000*00 

1 

.279*00 

2 

.2«2*Q0 

3 

.292*00 

4 

.316*00 

5 

, 36D+DU 

6 

. 3500*00 

1 

.263*00 

? 

.269*00 

3 

,289*00 

U 

. 32 7*00 

5 

. 344*00 

9 

. 4C,OD+ao 

1 

.351*00 

2 

.3*"8*J0 

3 

, 3 ? * n rj 

4 

• 472 + Of' 

5 

, 4 o[j + n: 

IG 

.4500*00 

1 

.459*00 

2 

.466*00 

3 

. 4 9 C ♦ C 

4 

.532*00 

5 

.590*00 

11 

. 5U00+G0 

1 

.541*00 

2 

.549*00 

3 

.574*00 

U 

.617*00 

5 

.679+00 

1 

.CLDO 

6 

. 100*01 

7 

. iro +01 

8 

. 100*01 

<3 

. irn*oi 

IQ 

, l-'Q + Ol 

2 

.5000-01 

6 

.9'^8 + OD 

7 

. 979*00 

8 

.946*00 

0 

. 96 1 *u^ 

1C 

. 9 P 5 • 0 0 

3 

. luPO+OQ 

6 

,644 *00 

7 

.656*00 

8 

.683+pn 

9 

, 9 3 1 + n f* 

ID 

. iro*Jl 

4 

.1500*00 

6 

. 726*00 

7 

. 7'"7*00 

e 

. 616*00 

0 

.91 o + or" 

ID 

. l'?4 + 01 

5 

. 2000*00 

6 

.599*00 

7 

. 0 ^ 6+0 0 

8 

. 75 3*nn 

0 

. 399*0'' 

10 

,1 10*01 

6 

.2500*00 

6 

.488*00 

7 

.575*00 

8 

, 709+cn 

9 

,399 + 0'' 

10 

.116*01 

7 

. 3000*00 

6 

.434*00 

7 

.544*00 

3 

,699 + 00 

9 

.91 3*00 

ID 

. 1->1*D1 

8 

. 3500+00 

6 

.466*00 

7 

. 576*00 

3 

. 7TC*0n 

9 

. 94 0*C0 

ID 

.124*01 

9 

.4000*00 

6 

.559*00 

7 

.661 *un 

e 

. 795*0"> 

9 

,979 + 1 '' 

10 

. 1 25*01 

lu 

.45C0+U0 

6 

.667+00 

7 

.762*00 

e 

, 0 79 ♦ 0 n 

0 

.10 3 + 01 

10 

.124*01 

11 

.6000*00 

6 

.7«^O + 0G 

7 

.659*00 

8 

. 9 73*0'" 

0 

. 109 + ul 

10 

,120*01 


i 

X position 

j 

p(i) 

j 

p(l) 

j 

p(l) 

3 

p(l) 

j 

p(l) 

1 

.GOOD 

1 

. l-'O + O 1 

2 

. l'"L*01 

3 

. i''0 + r 1 

4 

.100*01 

5 

. irc*oi 

2 

. 500G-01 

1 

.952*00 

2 

.951 +0T 

3 

.948*2'' 

4 

,945*00 

5 

.941+00 

i 

. iuoa*un 

1 

.657 + C C 

2 

.6'’5*JQ 

3 

. f 5 2 + r' 

4 

. 84 7*C^ 

5 

.842+UC 

4 

. 1500*00 

1 

.724*00 

7 

.773*20 

3 

.7 19*0'" 

4 

.714 + .'" 

5 

.710*00 

5 

. 20DD+00 

1 

.567*00 

2 

. 565*00 

3 

, 561 +2'" 

4 

•556+or 

5 

.552*00 

6 

.2500*00 

1 

,398*t r 

2 

.3^6* JO 

5 

, ' c 1 + r 

4 

. 3=4*00 

5 

. 378+UO 

7 

, 3oC0*00 

1 

.231*00 

2 

.229*00 

3 

. <^7l*^'' 

4 

.21 1+00 

5 

.2"0*DC 

6 

, 3S3C+C0 

1 

,613-01 

2 

. 774 - 0 1 

3 

.664-0] 

4 

,494-01 

5 

.289-01 

9 

.4000*00 

1 

-.401-01 

7 

-.4=7-01 

3 

-.67!'-'^! 

4 

- .082-0] 

5 

123*''b 

10 

.4500*00 

1 

-. 124*00 

2 

-.171*00 

3 

-.154+1''' 

4 

191 

5 

-.242*00 

11 

. 5000*00 

1 

-.164+OC. 

2 

-.174+00 

3 

-.2n2*(V" 

4 

-.251*'J'" 

5 

-.319*00 

1 

.0000 

6 

,i"’n*ci 

7 

. ira*ai 

8 

. ITQ.ni 

0 

.lOC + Gl 

IG 

. l'”U*Ql 

z 

. SoDO-Ol 

6 

.9 78*00 

7 

.937*00 

8 

, 940 + D'" 

9 

.95 1* jO 

10 

.968*00 

3 

. lonc+LO 

6 

.6 79 + CU 

7 

.643*00 

3 

.c. = 6*C'' 

9 

. « 8 5 * 2 0 

10 

,9 73* U’ 

4 

. is-ic+oo 

6 

,710*00 

7 

.719*00 

8 

, 7a4 ♦C'' 

9 

,794*oO 

1C 

.8=1*00 

5 

. 2000+00 

6 

,5'=4*CG 

7 

.567+00 

8 

. b r- r ♦ n 

9 

. 666*'j0 

10 

. 744+30 

b 

.2500*00 

6 

. 37S + 0C 

7 

.3«0*00 

0 

, 4 •> 2 ♦ f' '" 

0 

, a p 9 * f I' 

1-3 

.614*''u 

7 

. 3Lno*on 

6 

.193*00 

7 

.194*00 

a 

1 s+rf 

Q 

,26 0*2'' 

10 

.356*rc 

6 

. 3500*00 

6 

.799-02 

7 

-.897-0? 

5 

i7r-Di 

9 

- .99fc-t? 

10 

.201-31 

9 

.4^00*00 

6 

-.163*00 

7 

-.297*30 

8 

-. 6 5 1 ♦n'l 

9 

- . 297*1 r 

1C 

-.3=3*00 

10 

.4500*00 

b 

3''7*00 

7 

-.3«3*30 

8 

-.470+00 

C 

- .57 1 ♦O'" 

ir 

-.7'"3*0U 

11 

. SoDO+OO 

6 

-.430*00 

7 

523+on 

e 

-.6 '■7*‘'0 

0 

- , 439*, 0 

1C 

-.964*21 
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i 

X position 

j 

p(2) 

i 

p(2) 

j 

p(2) 

j 

p(2) 

i 

p(2) 

1 

• QQOO 

1 

.QOO 

2 

• one 

3 

.coo 

4 

.cno 

5 

.□ro 

2 

.5u00-01 

1 

.221-01 

2 

.198-01 

3 

.128-ni 

4 

.906-03 

5 

-.164-01 

3 

•1000400 

1 

.301-01 

2 

.256-01 

3 

.118-01 

a 

-.117-01 

5 

-.4‘'7-21 

4 

. 1500+00 

1 

.170-01 

2 

.114-01 

3 

-.8"9-0? 

4 

-.414-Jl 

5 

-.397-01 

B 

•2000+00 

1 

-.179-01 

2 

-.256-01 

3 

-.492-01 

4 

-.697-01 

5 

-.148+00 

6 

•2S00+00 

1 

-.770-01 

2 

-.854-01 

3 

-. 111 +cn 

4 

-.155 + Cjr 

5 

-.2?C + i-?b 

7 

.3000+00 

1 

-.156+00 

2 

-.165+00 

3 

-.141+00 

U 

-.235 + rjn 

5 

-.3''C + 0C 

6 

.3500+00 

1 

-.250+00 

2 

-.250+00 

3 

-. 2 ® 2 +cn 

4 

-.323+-r 

5 

-. 3“3+CC 

9 

.4000+00 

1 

-.340+00 

2 

-.355+00 

3 

-.377+Cr 

4 

-.41 3 + .-^ 

5 

-, 4 6b + '.\ 

1C 

•4500+00 

1 

-.441+00 

2 

-.447+00 

3 

-,466+Dn 

4 

-.496+0'’ 

5 

-,538 + fU 

11 

•5000+00 

1 

-.515+00 

2 

-.571+00 

3 

-.577+nr 

4 

-.564+ )r 

5 

-.509 + C.. 

1 

• 0000 

6 

.000 

7 

*Ul?0 

8 

.too 

9 

.HOC 

iC 

.cru 

2 

.5000-01 

6 

-.395-01 

7 

-.605-01 

a 

-.103+00 

9 

-.142*0'’ 

10 

-.132+ri; 

3 

.1000+00 

8 

-.909-01 

7 

-.146+00 

8 

-.215+00 

9 

-.289 + 1(0 

10 

3*2+c. 

4 

. 1500+00 

6 

-.154+00 

7 

-.236+00 

8 

-.334+rn 

Q 

-.444+Or 

10 

-.5''7+rL 

5 

.2000+00 

6 

-.220+QC 

7 

-.329+00 

e 

-.455+r- 

a 

-.6n4+ jr 

10 

-.770+ t . 

6 

•2500+00 

6 

-.3C8+QC 

7 

-.423+00 

8 

-.570+0? 

9 

-.755+cr 

1C 

-.9F2 + C . 

7 

.3000+00 

6 

-.309+00 

7 

-.Sf’e+cn 

3 

-. 666+On 

9 

- .875+ jr 

ID 

-.1 16+ni 

8 

• 3500+00 

6 

-.466+00 

7 

-.578+00 

8 

-.7^0+00 

9 

- .94 n+ J'* 

10 

-,1?4 + 'j1 

9 

•4000+00 

6 

-.535+00 

7 

-.628+00 

8 

-.7'4+rn 

9 

-.932 + ijr 

10 

-,l?t. + lil 

10 

.4500+00 

6 

-.592+00 

7 

-.649+00 

8 

-.742 + 0'' 

9 

-.853+Or 

1C 

ir2+ui 

11 

.5000+00 

6 

-.640+00 

7 

-.602+00 

8 

-.717+Cn 

9 

-.735 + 0'’ 

10 

-.7?i+ni; 


i 

X position 

j 

u<« 

J 

u<« 

j 


j 

u<« 

J 

u'« 

1 

.0000 

1 

.257+01 

2 

.267+01 

3 

.296+01 

4 

.348+31 

5 

.479+01 

2 

.5000-01 

1 

.242+01 

2 

.250+01 

3 

.276+ni 

4 

.321+01 

5 

. 3''U + C1 

3 

. ICOO+QO 

1 

.223+01 

2 

.2’0+01 

3 

.249+0! 

4 

.242+ul 

5 

.372+01 

4 

. 1500+0C 

1 

.2n2+Cl 

2 

.2r5+Dl 

3 

.217+01 

4 

.236+31 

5 

.?f 3 + 'j 1 

5 

.2600+00 

1 

. 177 + Cl 

2 

.178+01 

3 

.142+ni 

4 

. 1 86 +u 1 

c 

. 1=C + 'J1 

6 

.2500+00 

1 

.149 + 01 

2 

.148+01 

3 

. 14S + ? 1 

4 

.137+31 

S 

. l?L+'’l 

7 

.3JC0+00 

1 

.118+01 

2 

.116+01 

3 

.i''7+ni 

4 

.S92 + ,,r 

5 

.591+l't 

6 

.3500+00 

1 

.831+00 

? 

.7C2+03 

3 

.67C + m 

4 

.443+ jO 

5 

.775-01 

9 

.4000+00 

1 

.420+0C 

2 

.377+00 

3 

.245 + 0'’ 

4 

.592-?? 

5 

-.370 + '.', 

10 

.4500+00 

1 

-.617-01 

2 

-.l ''2 + 00 

3 

-.2?7+pn 

4 

- .447+’r 

5 

-.776 + 0 ., 

11 

.5000+00 

1 

-.618 + 00, 

2 

-.6‘=4+00 

3 

-. 7'^2 + f'’ 

4 

-.944+ jr 

5 

-.170+01 

1 

.O’uno 

6 

.548+01 

7 

.7'’i+ni 

8 

.970 + 0 1 

9 

.173+0? 

IL 

.145+:? 

2 

.5000-01 

6 

.489+01 

7 

.629+01 

8 

.tt?5+Cl 

9 

. nr + 3? 

in 

.14^+22 

5 

. lOOO + CO 

6 

.4PC+01 

7 

.4°2+0l 

8 

.613+01 

P 

.773+31 

Iv. 

.9-’6 + u1 

4 

.1500+00 

6 

.295+01 

7 

.332+31 

a 

. 369+DI 

9 

. 596+01 

in 

.4''C*ri 

5 

. 2000+00 

6 

.187+01 

7 

.172+01 

a 

.136 +C1 

9 

. 363+ jO 

1C 

-. 14Q+0 1 

6 

.2500+00 

6 

.894+00 

7 

.327+00 

8 

-.648+rn 

0 

- .244+01 

ID 

-.565+ni 

7 

. 300D+00 

6 

.923-01 

7 

-.722+00 

0 

-.2"'5 + 01 

9 

-.421 + fil 

10 

-.74 ., + 01 

8 

.3500+00 

6 

-.506+00 

7 

-.179+01 

p 

-.273+C 1 

9 

-.477 + 1.1 

iO 

-.7'’!, + 01 

9 

.4u00+C0 

6 

-.930+00 

7 

-.174+31 

8 

-.286+01 

0 

- .445+ul 

1C 

6 4 4 + n 1 

10 

.45ra+oo 

6 

-.125+01 

7 

-.149+01 

0 

-.272+01 

9 

-.374+Cl 

1C 

-.4 75 + '^ 1 

11 

.5UOC+00 

6 

-.155+01 

7 

-.i'>7 + Dl 

8 

-.<;4S*n] 

0 

3nc+ji 

10 

-.3'3 + -)l 


i 

X position 

j 

u<2) 

J 

„(2) 

j 

u<2) 

j 

u(2) 

j 

„(2) 

i 

.0000 

1 

-.6r2+uC 

2 

-.640+00 

3 

-,958+Cn 

4 

-.11 7+(il 

5 

-. 14V+01 

2 

,5uno-oi 

1 

-.191+01 

2 

-.197+01 

3 

-.213+01 

4 

- .244+ul 

5 

-.?91+L 1 

3 

,1600+00 

1 

-.2PS+C1 

2 

-.291+01 

3 

-.312+01 

4 

-.348+01 

5 

-.4' 6 + '“ 1 

4 

.1500+00 

1 

-.356+01 

2 

-.563+01 

3 

345+01 

4 

-.4?4+„] 

5 

- . 4 >• : + ‘ 1 

b 

. 2 coo+ac 

1 

-.400+01 

2 

-,4'’7+ai 

3 

-.429 + C 1 

4 

-.467+ 11 

5 

-.5?b+?l 

6 

.2500+00 

1 

-.418+01 

2 

424 + U 1 

3 

-, 444+01 

4 

-.478+ :i 

5 

-.5 3,. + 1.1 

7 

.3000+00 

1 

-.408+01 

2 

-.414+01 

3 

-.431+01 

4 

-.461*31 

5 

-.5''4+ni 

8 

.3500+00 

1 

-.376+01 

2 

-.541+31 

3 

-, 395 +ri 

4 

-.*t?r’ + oi 

5 

-.4^3+.'i 

9 

. 4L.OO+CO 

1 

-.325+01 

2 

-.329+01 

3 

-.341+01 

4 

-.360+01 

5 

746+' 1 

10 

.4500+00 

1 

-.262+01 

? 

-.266+01 

3 

-.2 75 + Cl 

4 

- .289+LI 

5 

i'-s + O 1 

11 

.5000+00 

1 

-.194+01 

2 

-.106+01 

3 

-.2"2+ri 

4 

-.212 + i:l 


?'‘6+ ■ 1 

1 

.CLOU 

6 

-.196+01 

7 

-.265+01 

8 

566+ri 

9 

- . 5 1 6 + u 1 

10 

-.74« + : 1 

2 

.5000-01 

6 

-,363+ni 

7 

-.47[ +Q1 

8 

-.b^o + p 1 

9 

-.872+Jl 

12 

- .124 + ' 2 

3 

.1000+00 

6 

-.403 + C 1 

7 

-.672+Cl 

8 

-.618+Cl 

Q 

-.112+U? 

lO 

- . 1 ' a ♦ 0 2 

4 

.1500+00 

6 

-.576+01 

7 

-.712+01 

8 

vl4+?l 

0 

-.17?+0? 

1'’ 

- . 1 7 J + f. ? 

5 

.200D+D0 

6 

-.612+01 

7 

-.?r6*oi 

a 

-.91 fc+Pl 

C 

- . UP + u7 

in 

-.14-,+ '2 

6 

.2500+00 

6 

-,6''3+tn 

7 

-.7'’3+Cl 

« 

-.839+01 

9 

- .in?+u? 

1C 

I’n + r-? 

7 

.3O00+CO 

6 

-.560+01 

7 

-.670+01 

a 

71 3*ri 

P 

-.302+01 

1C 

-.rtPP + H 

a 

. 3500+00 

6 

-. 493 +Cl 

7 

-.556+01 

s 

-.572+01 

0 

-.58 1+i: 1 

10 

-.5-1+01 

9 

.4000+00 

6 

-.414+01 

7 

-.439+01 

8 

-,447 + n 

9 

-.41 1+ul 

1C 

-.276+L 1 

10 

,4500+CO 

6 

-.329+01 

7 

-.347+01 

8 

-.349+01 

0 

-.316+Jl 

1C 

- • 2 ■' 0 + •'■ 1 

11 

.5000+00 

6 

-.241+01 

7 

-.2‘^7 + 01 

8 

-.27C+n 

0 

-.277+Cl 

1C 

-.27? + {’ 1 


i 

X. position 

i 

v<« 

j 

yd) 

j 

yd) 

i 

yd) 

j 

yd) 

1 

.0000 

1 

-.159-02 

2 

.i«D+on 

3 

.383 + C'’ 

4 

.621+rn 

5 

.9?6 + ni' 

2 

.5000-01 

1 

-.239-02 

2 

.340+00 

3 

.721*m 

4 

.117+51 

5 

. I7b+ni 

3 

. 1000+00 

1 

-.276-02 

2 

.472+00 

3 

.997* rn 

4 

.161+Jl 

5 

.2’9 + fjl 

4 

. 1500*00 

1 

-.264-02 

2 

.564+00 

3 

. iie+n 

4 

.I9C + 1, I 

5 

.279+n 

5 

.2000+00 

1 

-.205-02 

2 

.611+00 

3 

.l'’7*01 

4 

.292+^1 

5 

.293+01 

6 

.2500+00 

1 

-.115-02 

2 

.614+00 

3 

.127*01 

4 

.199+LI 

5 

.2ft3 + C 1 

7 

.3L00+C0 

1 

-.170-03 

2 

.500+CC 

3 

. 118*01 

4 

.183*01 

5 

.296+C 1 

6 

.3500+00 

1 

.662-03 

2 

.520+00 

3 

,inB*r 1 

4 

.159.C1 

b 

.217*ill 

9 

.4U0D+00 

1 

.120-02 

2 

.446*00 

3 

,e«7*r'’ 

4 

. 133 + Ll 

5 

. 1 75+01 

10 

,4500+00 

1 

.145-02 

2 

.371+0? 

3 

.7‘?6*P'' 

4 

.n6+ci 

5 

.1-^6 + Cl 

11 

.5000+00 

1 

.149-02 

2 

.3'-2*30 

3 

.585*0'’ 

4 

.842 + J*" 

5 

. i''4+ni 

1 

.0000 

6 

.135*01 

7 

.197+01 

5 

.291*01 

Q 

.435+Cl 

10 

.b64+01 

2 

.5000-01 

6 

.254*01 

7 

.366+01 

8 

.531+ri 

0 

. 7 ft (, + L 1 

1C 

.116+12 

3 

.1000+00 

6 

.343*01 

7 

.489+01 

8 

.699. PI 

9 

.101+02 

k; 

. 149+C 2 

4 

.1500*00 

6 

.395*01 

7 

.553+01 

8 

.777*01 

9 

.iin+0? 

10 

. 16U + 02 

5 

.2000*00 

6 

.407*01 

7 

.558*01 

8 

.763*01 

9 

.105*02 

10 

. 1“7*C2 

6 

.2500*00 

6 

.385*01 

7 

.511+01 

B 

.672*01 

9 

.803+Cl 

10 

.117*02 

7 

. 3000+00 

6 

.337*01 

7 

.430+01 

8 

.535+01 

9 

.650+01 

ID 

.767 + 01 

8 

.3500*00 

6 

.275*01 

7 

.333*01 

B 

.384*01 

9 

.412+01 

10 

.381+01 

9 

.4000*00 

6 

.212*01 

7 

.240*01 

8 

.248+01 

9 

.216*01 

10 

.985*00 

lU 

.4500*00 

6 

.157+01 

7 

.163*01 

6 

. 145+ni 

0 

.857*00 

10 

-.4U5+3C 

11 

.5000*00 

6 

.115+01 

7 

.109*01 

8 

.735*00 

0 

.129+30 

10 

-.981+00 


50 



i 

X position 

3^ 

y(2) 
V2 a 

3 


3 

V® 

3 


3 

v(2) 

''a 

1 

.0000 

1 

-.235-02 

2 

.369*00 

3 

.790*00 

4 

.130*01 

5 

.197*01 

2 

.5000-01 

1 

-.218-02 

2 

.355*00 

3 

.761*00 

4 

.126*01 

5 

. 190*01 

3 

. 1000*00 

1 

-.182-02 

2 

.342*00 

3 

.731*00 

4 

.120*01 

5 

.182*01 

4 

.1500*00 

1 

-.128-02 

2 

.329*00 

3 

.699*00 

4 

.114*01 

5 

.171*01 

5 

.2000*00 

1 

-.622-03 

2 

.313*00 

3 

.658*00 

4 

.106*01 

5 

.156*01 

6 

.2500*00 

1 

.496-04 

2 

.292*00 

3 

.606*00 

4 

.963*00 

5 

.139*01 

7 

.3000*00 

1 

.619-03 

2 

.266*00 

3 

.542*00 

4 

.845*00 

5 

.118*01 

8 

. 3500*00 

1 

.101-02 

2 

.236*00 

3 

.473*00 

4 

.719*00 

5 

.969*00 

9 

.4000*00 

1 

.124-02 

2 

.208*00 

3 

.410*00 

4 

.605*00 

5 

.779*00 

10 

.4500*00 

1 

.138-02 

2 

.190*00 

3 

.367*00 

4 

.526*00 

5 

.646*00 

11 

.5000*00 

1 

.145-02 

2 

.183*00 

3 

.352*00 

4 

.498*00 

5 

.598*00 

1 

.0000 

6 

.289*01 

7 

.419*01 

a 

.607*01 

9 

.882*01 

10 

.129*02 

2 

.5000-01 

6 

.279*01 

7 

.405+01 

8 

.587*01 

9 

.853*01 

10 

.126*02 

3 

.1000*00 

6 

.265*01 

7 

.381*01 

8 

.548*01 

9 

.792*01 

10 

.116*02 

4 

.1500*00 

6 

.245*01 

7 

.347*01 

8 

.489*01 

9 

.692*01 

10 

.996*01 

5 

.2000*00 

6 

.220*01 

7 

.303*01 

8 

.411*01 

9 

.557*01 

10 

.759*01 

& 

.2500*00 

6 

.189*01 

7 

.250*01 

8 

.321*01 

9 

.400*01 

10 

.484*0 1 

7 

.3000*00 

6 

.155*01 

7 

.193*01 

8 

.227*01 

9 

.243*01 

10 

.213*01 

8 

.3500*00 

6 

.120*01 

7 

.138*01 

8 

.141*01 

9 

.108*01 

10 

-.693-01 

9 

.4000*00 

6 

.901*00 

7 

.918*00 

8 

.720*00 

9 

.978-01 

10 

-.142*01 

10 

.4500*00 

6 

.689*00 

7 

.597*00 

8 

.270*00 

9 

-.450*00 

10 

-.187*01 

11 

.5000*00 

6 

.611*00 

7 

.477*00 

8 

.106*00 

9 

-.614*00 

10 

-.182*01 


Problem 4 


Plane 



1 = 0 {use e;^ = oo} 
(d) 


77 = 0. 6 
7?+ = 0. 6 
M = 0. 5 

= 0.071 - 0. 151 (upper) 
= IXIO^® (centerline) 
P(0, y) = 1 


n = 10 
m = 10 


51 


The calculated values are 


Ax = 0. 5 

Ay = 0. 05556 

Eq = 2. 18 

= 0. 853 (x = 0. 5) 
AdB = -4. 07234 


k 








Ok 


1 

-.926*00 

-.233*00 

-.100*01 

.000 

.385*01 

.ncD 

-.100*01 

.COC 

-.926+00 

2 

.132*02 

-.118*02 

-.200+01 

.000 

-.291*02 

. 1S5*'’2 

.QOD 

.000 

.169+02 

3 

-.926*00 

-.233*00 

.000 

.000 

.385+01 

.304-29 

-.200+01 

.OOC 

-,92t + r,0 

•* 

-.185*01 

-.465*00 

-.ICO+Cl 

-.796*01 

.385+01 

. 164 + "2 

-.IDU+Ol 

-.796*01 

.COP 

5 

.282*02 

-.195*02 

-.200+01 

-.159*u2 

-.310+02 

.279+02 

.uno 

.cno 

.POD 

6 

-. 185*01 

-.465*00 

.coc 

.uon 

. 385+01 

.164+02 

-.200+01 

-.159+02 

.nor 


i 

X position 

i 

|P| 

j 

|Pl 

j 

|Pl 

i 

|P| 

j 

|Pl 

1 

.0000 

1 

. ino*oi 

2 

.ino+Di 

3 

.ino+Pi 

4 

,100+01 

5 

.100*01 

2 

. SuOO-Ol 

1 

.942*00 

2 

.941+00 

3 

. 940+rn 

4 

.937+00 

5 

.935+00 

3 

. lDDO+00 

1 

.843+00 

2 

.843+00 

3 

.640+00 

4 

.838+QP 

5 

.8'^6 + nc 

4 

. 1500+00 

1 

.715+00 

2 

.714+00 

3 

.713+00 

4 

.712+UC 

5 

,716+OC 

5 

.2000+00 

1 

.572+00 

2 

.572+00 

3 

.573+00 

4 

.578+ur 

5 

.59Q+DC 

6 

.2500+00 

1 

.441*00 

2 

,443+On 

3 

.448*0'' 

4 

.461+OP 

5 

,4P7+0D 

7 

. 300D+D0 

1 

.366+00 

? 

.37u*on 

3 

.3"2+rr 

4 

.405+0r 

5 

, 4*»4 + 0 u 

a 

.35DD+O0 

1 

.384+00 

2 

.389+00 

3 

.406+0" 

4 

.436+00 

5 

.482+00 

9 

.4o0C+00 

1 

.466*00 

2 

.473+00 

3 

,492+00 

4 

.526*00 

5 

,574+00 

lU 

.4500+00 

1 

.560+00 

2 

.567*00 

3 

.589+no 

4 

.625+00 

5 

.676+00 

11 

.50DD+U0 

1 

. 632 + 00 

2 

.64Q+Da 

3 

.663+Cn 

4 

.7n2+uC 

5 

,753+00 

1 

.GOOD 

6 

.inO+01 

7 

.inc+Dl 

8 

.100+01 

9 

. 100 + 01 

ID 

. IDO + Dl 

2 

.5000-01 

6 

.932+00 

7 

.932*00 

8 

.933+PO 

9 

,936+00 

10 

,9ua+00 

3 

. luOD*CO 

6 

.8’B+OO 

7 

.845+00 

8 

,860+00 

Q 

.833*jr 

10 

.910+00 

4 

. 15P0 + D0 

6 

.727+00 

7 

.790+00 

3 

.790+00 

O 

,348*00 

10 

.916*00 

5 

.2000*00 

6 

.616*00 

7 

.662+00 

3 

.734+CC 

9 

.833+00 

ID 

,9S6+nC 

6 

.2500*00 

6 

,531+00 

7 

.601+00 

8 

.7"3+PO 

c 

,841+On 

10 

.102+01 

7 

. 3000*00 

6 

,504+00 

7 

.689+00 

8 

.7P8+nn 

9 

,870+00 

10 

,1P9+C1 

8 

.3500*00 

6 

.547+00 

7 

.634+00 

8 

.753+00 

9 

.918+00 

ID 

. 1 16+01 

9 

,4000+00 

6 

.638+00 

7 

.771+00 

0 

.63D+P" 

o 

.981+00 

10 

.171*01 

10 

.4500*00 

6 

.743+00 

7 

.825+CO 

8 

.927+PO 

9 

.106+01 

10 

,125+01 

11 

.5000+00 

6 

.832+00 

7 

.974+00 

8 

.IP3+PI 

9 

.115+01 

10 

.125+01 


.233*DO 

-.527+01 

.233*00 

.000 

.000 

.000 


i 

X position 

i 

p(« 

3 

p(i) 

]■ 

p(i) 

3 

P<« 

3 

p(l) 

1 

.0000 

1 

.100+01 

7 

,1^0+01 

3 

.100+01 

4 

.inc+Dl 

5 

.IPQ+Ol 

2 

.5000-01 

1 

,934+DD 

2 

.934*00 

3 

.933+00 

4 

.931+un 

5 

.93D+DU 

3 

. IDOO+CO 

1 

.e?D+oc 

2 

.870+00 

3 

.8?0 + 0<^ 

4 

.8?i*'jr 

5 

.824+00 

4 

.isno+up 

1 

,630*00 

2 

.6«1*DC 

3 

, 6 B 4 ♦ m 

4 

.69P+or 
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TABLE I. - COEFFiaENTS IN DIFFERENCE EQUATION (EQ. (54)) 
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TABLE n. - SUMMARY OF OPTIMIZATION 
RESULTS (REF. 9) 


[Length-height ratio, L*/H*, 3. 43; Mach 
number, M, 0. 3; plane-wave source.] 


Dimensionless 

frequency, 

V 

Specific 

acoustic 

impedance, 

? 

Sound 

attenuation, 
A dB, 
dB 

0. 5 

0. 31 - 0. 1 i 

70.6 

1.0 

.65- .5 1 

39. 2 

1.535 

.85 - 1.0 i 

22.6 

2.0 

. 90 - 1. 4 i 

14.3 

5.0 

. 85 - 2. 45 i 

4.0 

10.0 

1. 40 - 3. 6 i 

2. 30 



Dimensionless axial coordinate, x 


Figure 1. - Typical pressure profile for sound propagation in 
two-dimensional, soft-wall duct. Dimensionless frequency, 
ri, 1; dimensionless duct length, L*/H* 3. 



Dimensionless axial coordinate, x 


Figure 2. - Sensitivity of real pressure component p^^^ to choice of wave enve- 
lope of frequency for two-dimensional hard-wall duct. Dimensionless 
frequency, r\, 1; dimensionless duct length, L"/H'" = 3. 
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Figure 3. - Effect of axial duct length and dimensionless frequency 
on attenuation at optimum impedance in two-dimensional, soft- 
wall duct assuming that wave envelope frequency equals t?. 
(From ref. 4. ) 



grid points 

Figure 4. - Coordinate and grid-point representation of two-dimensional, soft-wall duct. 





. A.A JA fV-l . 


Jl) 


(a) Conventiona! pressure solution. 

p-i^px 


^-A-A-A-A-X _ 
\i\JALU^Ui 


2 

Dimensionless axial coordinate, x 


(b) Wave envelope solution. 

Figure 6. - Typical pressure profiles for sound propagation in two- 
dimensional. soft-wall duct for dimensionless frequency ri of 1. 
dimensionless duct length L"/H* of 3, and Mach number M 
of 0. 5. 
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Axial position along upper wall. 

(a) Exit condition of reference 3. Sound attenuation. (b) Exit condition of this report. Sound attenuation. 
AdB, 6.1 decibels. AdB, 3.4 decibels. 

Figure 7. - Comparison of dimensionless transverse acoustic particle velocity v at wall gs calculated from momen- 
tum equation and impedance condition. Dimensionless frequency, rj, 0.6; dimensionless duct length, L"/H* 
0.5i Mach number, M, 0.5; specific acoustic impedance, C 0.071 - i 0.151. 



Figure 8. - Effect of dimensionless frequency on attenuation at 
optimum impedance in two-dimensional soft -wall duct. Dim- 
ensionless duct length, 3.43; Mach number, M, 0.3. 
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